SPIN, TWIST AND HADRON STRUCTURE 
IN DEEP INELASTIC PROCESSES* 



R. L. Jaffe 

Center for Theoretical Physics 
Massachusetts Institute of Technology 
Cambridge, MA 02139 U.S.A. 

and 

Department of Physics 
Harvard University 
Cambridge, MA 02138 U.S.A. 



Notes for Lectures Presented at the 
International School of Nucleon Structure 
The Spin Structure of the Nucleon 
Erice, 3-10 August 1995 



ABSTRACT 

These notes provide an introduction to polarization effects in deep inelastic pro- 
cesses in QCD. We emphasize recent work on transverse asymmetries, subdominant 
effects, and the role of polarization in fragmentation and in purely hadronic processes. 
After a review of kinematics and some basic tools of short distance analysis, we study 
the twist, helicity, chirality and transversity dependence of a variety of high energy 
processes sensitive to the quark and gluon substructure of hadrons. 



*This work is supported in part by funds provided by National Science Foundation (N.S.F.) grant 
#PHY 92-18167 and by the U.S. Department of Energy (D.O.E.) under contracts #DF-FC02-94ER40818 
and #DF-FG02-92ER40702. 

MIT-CTP-2506 and HUTP-96/A003 January 1996 



Contents 
Introduction 

1 Kinematics and other Generalities 

1.1 Deep Inelastic Scattering 

1.1.1 Basic Variables 

1.1.2 Cross Section and Structure Function^ 

1.1.3 Structure li'imctioni 

1.1.4 Scaling and Kinematic Domains 

1.1.5 Flavor Ceneralizations 

1.1.6 Cross Section for Electron-Hadron Scattering 

1.2 Other Basic Deep Inelastic Processes 

1.2.1 Inchisive e"^ e~ Annihilation 

1.2.2 Inclusive e"^e~ Annihilation with One Observed Hadron 

1.2.3 Lepton Pair Production 



2 


Deep Inelastic Processes from a Coordinate Space Viewpoint] 


2.1 


e+e^ hadrons - The Short-distance Expansion 




1.2 


Ip ^ IX - 'I'he Light-Cone Expansion 






2.2.1 Light-Cone Coordinates and Formulation of Deep Inelastic Scattering 






2.2.2 Deep Inelastic Scattering and the Short Distance Expansion 




2.3 


e^e^ — > h X- Once Again, the Light-Cone 




2.4 


PP ^ l+l-X - The Drell-Yan Process 




2.5 


Dominant and Subdominant Diagrams! 



3 Deep Inelastic Scattering and Generalized Distribution Functions I 

3.1 Twist 

p. 1.1 Twist and the OPE| 



p. 1.2 Examples and a Working Redefinition of Twist 


B.1.3 Spin and Twist 


3.2 Dominant Diagram in Coordinate Space 




3.3 Learning from Light-Cone Quantization 


B.3.1 Light-Cone Kinematics 


p. 3. 2 Dirac Algebra on the Light-Cone 




p. 3. 3 Independent Degrees of Freedom 




3.4 i'he Parton Model 




p. 4.1 Twist-Two 




P.4.2 Twist-Thre( 


3 



4 


Deep inelastic Scattering and Generalized Distribution iYinctions 11 




4.1 Helicity Amplitudes 




4.2 Quark Distributions 


m 


'I'argets with Spin-U, 1/2 and 1 




4.2.1 Spin-0 Target 




12.2 Spin-1/ 2 Target 





1 



4.2.3 Spin-1 Targetl 43 



4.3 Quark Distribution .b unctions for the JNucleon| 43 



4.3.1 Nucleon Spin Structure at Twist-Two| 44 

4.4 Transverse Spin in QCD| 47 

4.5 Twist-Three: Physics with 92(3^; 49 

4.5.1 Operator Product Expansion] 49 

4.5.2 Wandzura - Wilczek Decomposition of .92 51 

4.5.3 I'he Burkhardt Cottingham Sum Rule] 53 



4.5.4 The Evolution of ^2! 54 

4.5.5 'I'he Shape ot" (72I 55 



5 The Drell-Yan Process 


5.1 Operator Analysis 




5.2 Polarized Drell-Yan: a Briet Summary 





6 Annihilation and Quark Fragmentation 



Functions 



56 



62 



62 



|6.1 Single Particle Inclusive Annihilation] 63 

3.2 Polarized q ^ A Frap;mcutatiou Functions from e'^e" — » A + X 



66 



|6.3 Observing hi[x,Q^) in Electroproduction] 67 



Introduction 

In recent years hadron spin physics has emerged as one of the most dynamic areas of 
particle physics. During the same period the field has got considerably more complicated. 
In times past only longitudinal asymmetries, that have simple parton model interpretations, 
attracted much attention; only dominant effects, that scale in the Bjorken limit, were experi- 
mentally accessible; and only relatively crude experimental data were available. Now interest 
has spread to transverse polarization asymmetries, subdominant effects, polarization effects 
in fragmentation and in purely hadronic processes. The aim of these lectures is to present an 
introduction to spin dependent effects at dominant and subdominant order in deep inelas- 
tic processes including deep inelastic scattering of leptons, e~^e~ annihilation, and Drell-Yan 
processes. The methods can be extended relatively straightforwardly to other spin dependent 
effects in hard processes. 

In a short set of lectures some detail and background must be sacrificed. As for back- 
ground, I will assume that readers are familiar with the elementary parton model treatment 
of highly inelastic processes in the "infinite momentum frame" . Anyone who is not familiar 
with basic parton model ideas should consult standard textbook presentations.&H Although 
I will have a lot to say about the parton model, it may look poorly motivated to someone 
who has not seen the ideas presented in their simplest form first. As for detail, I will mostly 
ignore the complications of QCD radiative corrections, normally included via the renormal- 
ization group. There are many excellent treatments including books by Collins,!! Mutai and 
most recently in a context particularly well suited to these lectures, by Roberts.! Of course 
radiative corrections and the momentum scale dependence they generate are central to the 
understanding of QCD. Some important aspects are covered in Al Mueller's lectures in this 
volume. Here we will be interested in the classification of scattering amplitudes in terms of 
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helicity, chirality, twist, etc. - a classification which is largely (but not entirely) independent 
of radiative corrections. In many cases the soft, InQ^ dependence they generate can be re- 
garded as decorations of our primary results. Where this is not the case, I will try to warn 
the reader and refer to the appropriate literature. 

The main question to be addressed here is: How can one classify and interpret the wide 
variety of spin dependent phenomena expected in hard processes? Which phenomena are 
displayed in which experiments? What are the selection rules enforced by the symmetries 
of QCD? Which phenomena dominate at large-Q^, which are suppressed, and what is the 
physical origin of the suppression? In short, the object is to provide the background for 
both experimental and theoretical analysis of spin effects in hard processes. In contrast, I 
will resist almost entirely the temptation to speculate about the origins of spin effects based 
on models of hadron structure. These notes are not intended to be an introduction to the 
so-called "spin crisis" which grew out of the observation that only a small fraction of the 
nucleon's spin is carried by the spin of quarks. Theorists will not find their own or my own 
favorite explanation of the spin crisis in these lectures. That is a subject for another school. 

Certain predictions of perturbative QCD are admired for being very general and inde- 
pendent of the difficult details of hadron structure. Examples include the cross section for 
e~^e~ — > hadrons, event shapes in e~^e~ annihilation, the InQ^ dependence of deep inelas- 
tic structure functions, and the Gross-Llewellyn Smith and Bjorken Sum Rules. Studies of 
these processes provide essential tests of QCD. These will not be major topics here. I will 
assume that perturbative QCD is correct and use it as a sophisticated probe of the poorly 
understood dynamics of confinement. As we shall see, perturbative QCD is by now so well 
understood that it is possible to "tune" the probe to measure the nucleon expectation values 
of a variety of quark and gluon distributions and correlations within hadrons. Probes can 
be selected for spin, twist and flavor quantum numbers, and can be used either to analyze 
the structure of hadronic targets or reaction fragments. No other approach yields such well 
defined information about hadronic bound states. This information may help guide us to a 
better understanding of confinement from first principles. 

Many aspects of these lectures are based on work performed in collaboration with 
Xiangdong Ji. The reader who wishes to explore subjects in greater depth should look at 
refs.,llrtj as well as other references provided in the text. I would like to thank Xiangdong 
for the pleasure of this long collaboration. Thanks are also due to Matthias Burkardt, Gary 
Goldstein and Aneesh Manohar who collaborated on other projects related to this work. 
In addition I have benefited greatly from discussions with Guido Altarelli, Xavier Artru, 
Ian Balitsky, Vladimir Braun, Gerry Bunce, John Collins, Vernon Hughes, Gerd Mallot, Al 
Mueller, Richard Milner, John Ralston, Phil Ratcliffe, Klaus Rith, Jacques Soffer, and Linda 
Stuart. 

These lectures grew out of talks at schools and conferences in the early 1990's. A version 
presented at the 1992 GraduiertenkoUeg of the Universities of Erlangen and Regensberg at 
Kloster Banz, Germany, was recorded by Drs. H. Meyer and G. Piller. The present version 
is based on a manuscript prepared by Drs. Meyer and Piller from their notes. I would like 
to thank them for the substantial work they undertook at that time. Subsequently I have 
edited, reformulated and expanded the notes, most recently for the Erice School on the 
Internal Spin Structure of the Nucleon. 
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1. Kinematics and other Generalities 



The organizers of the school asked if 1 would briefly introduce the kinematic and dy- 
namical variables common in the study of deep inelastic processes. So before getting down 
to the business of dynamics here is a short summary — the cogniscenti will certainly want to 
skip this section. Others, who may be familiar with less streamlined notation might wish at 
least to look at eqs. ( [Ll4D , (IlTsD , (11211) , and (\L2^) . I hope students with less background 



in perturbative QCD will find this section useful. 
1.1. Deep Inelastic Scattering 

1.1.1. Basic Variables 

Deep inelastic scattering (DIS) is the archetype for hard processes in QCD: a lepton — 
in practice an electron, muon or neutrino — with high energy scatters off a target hadron 
— in practice a nucleon or nucleus, or perhaps a photon — transferring large quantities 
of both energy and invariant squared-four-momentum. For charged leptons the dominant 
reaction mechanism is electromagnetism and one photon exchange is a good approximation. 
For neutrinos either (charged current) or (neutral current) exchange may occur. The 
weak interactions of electrons may also be studied either by means of small parity violating 
asymmetries originating in 7 — Z*^ interference, or by means of the charged current reaction 

e" Ue. 

We are primarily interested in experiments performed with polarized targets. Neutrino 
scattering experiments require far too massive targets for polarization to be a practical 
option, so we will ignore them, although l^-exchange has been observed in e~p + X, 
at HERA,lli and could be extended to a polarized target, at least in principle. Thus we 
are mainly limited to charged lepton scattering by one photon exchange. The kinematics is 
shown in fig. (|^). The initial lepton with momentum k and energy E exchanges a photon 
of momentum q with a the target with momentum P. Only the outgoing electron with 
momentum k' and energy E' is detected. One can define the two invariants 

q^ = {k- k'f = ql-q^ = -AEE' sm\ei2) = -Q^ < q (1.1) 
v = P-q = M{E-E'), (1.2) 




Fig. 1. Kinematics of lepton-hadron scattering in the target rest frame 



4 



where the lepton mass has been neglected (and will be neglected henceforth). The meaning of 
the scattering angle 6 is clear from fig. (Q). Unless otherwise noted, E', 9 and = E — E' 
refer to the target rest frame. The deep inelastic, or Bjorken limit is where and u both 
go to infinity with the ratio, x = Q"^ /2v fixed, x is known as the Bjorken (scaling) variable. 

Since the invariant mass of the hadronic final state is larger than or equal to the mass 
of the target, (P + g)^ > M^, one has < x < 1. It is convenient also to measure the energy 
loss using a dimensionless variable, 

0<y=-^<l. (1.3) 

- ^ ME - ^ ' 

We will find E, Q^, x, and ?/ to be a useful set of variables. Note that it is overcomplete 
since xy = Q'^/2ME, and note also that what we define as u differs from common usage by a 
factor of M. The behavior of cross sections at large is much more transparent using these 
variables than using the set {E,E', 9) favored by experimenters for the reason that ^ as 
— > oo at fixed x and y. 

1.1.2. Cross Section and Structure Functions 

The differential cross-section for inclusive scattering (eP —>■ e'X) is given by: 

da = - \, yT] L ,f \AW2'nf6\P + q-yp,). (1.4) 

The flux factor for the incoming nucleon and electron is denoted by J = 4P-k, which is equal 
to J = AME in the rest frame of the nucleon. The sum runs over all hadronic final states 
X which are not observed. Each hadronic final state consists of nx particles with momenta 
Pi (I]"=iPj = px)- The squared-amplitude can be separated into a leptonic {l^''') and a 
hadronic {W'^'^) tensor (see fig. @): 

A 2 «2 



An 



- Q-J^'W.u, (1.5) 

where a ~ 1/137 is the electromagnetic fine structure constant. The leptonic tensor l^'^ is 
given by the square of the elementary spin 1/2 current (summed over final spins): 

= J2u{k,s)Yu{k',s')u{k',s')Yu{k,s) 

s' 

= 2{k'''k'' + k'-'k") - 2gf"'k ■ k' + 2ie'"'^''qxs^, 

(1.6) 

and consists of parts symmetric and antisymmetric in /i and u. The antisymmetric part is 
linear in the spin vector s, which is normalized to = — m^. While the leptonic tensor is 
known completely, W'^'^ , which describes the internal structure of the nucleon, depends on 
non-perturbative strong interaction dynamics. It is expressed in terms of the current J'^ as: 

AttW^"" = Y.^PS\.P\X){X\.r\PS){2TxY5{P + q~px) (1.7) 



X 
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Fig. 2. The squared amplitude A for electron-hadron scattering can be separated into a leptonic tensor 
l'^'^ and a hadronic tensor W'^'^. 



The steps leading from eq. ( |1.7| ) to eq. ( [L.8D include writing the 5 function as an exponential, 

{27r)*6\K) = J dV^-^, (1.9) 

translating the current, e'^-^^-P^\PS\Jf'{0)\X) = {PS\J^'{^)\X) , and using completeness, 
Y,x = 1- Note that another term has been subtracted to convert the current product 

into a commutator. It is easy to check that the new term vanishes for g*^ > which is the case 
for physical lepton scattering from a stable target. The subscript c means that the graphs 
associated with the matrix element must be connected. Finally, note that the states are 
covariantly normalized to: 

{P\P') = 2E{2nf5^{P - P'). (1.10) 



The optical theorem: 



2txW^'' = Im T'"' 



;i.ii) 



relates the hadronic tensor to the imaginary part of the forward virtual Compton scattering 
amplitude, T: 

(1.12) 



T^. = W d^^e''^<{PS\T{J^{OM0))\PS) 



as shown graphically in fig. (|]). 



1.1.3. Structure Functions 

Using Lorentz covariance, gauge invariance, parity conservation in electromagnetism and 
standard discrete symmetries of the strong interactions, W^" can be parametrized in terms 
of four scalar dimensionless structure functions Fi{x, Q^), F2{x, Q^), gi{x, Q"^) and g2{x, Q"^). 
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Im 




Fig. 3. The optical theorem relates the hadronic structure tensor, W^"^ to the imaginary part of forward 
{P = P' , q = q'), virtual {Q^ < 0) Compton scattering 



They depend only on the two invariants and u, or ahernatively on and the dimension- 
less Bjorken variable x. Splitting W^" into symmetric and anti-symmetric parts we have, 



with 




VT'^'^ = H^W + (1.13) 

(1-14) 
(1.15) 

where 5"°" is the polarization vector of the nucleon (S"^ = — M^), P ■ S" = 0. S''^ is a 
pseudovector. Since W^^^^ is a normal tensor, parity demands that the appear with 
another pseudotensor, and the only one available is the e^^'^^ , Students often ask why 
W^^ depends only linearly on - what is wrong with S'^S'^, for example? Lorentz in- 
variance demands that W'^'^, defined in eq. ( p..8|) be linear in the initial and final nucleon 
spinors, U{P,S) and U{P,S). Tensors constructed from these are either spin independent 
{U{P, S)YUiP, S) = 2P^) or linear in 5^ ((f7(P, S)Yl5U{P, S) = 23"), but that is the end 
of it. 

Note also that W^" is dimensionless (we shall have more to say about operator di- 
mensions shortly). Factors of z/ have been judiciously introduced into eqs. (|1.15|) and ( p.. 141) 
so that the four structure functions. Pi, P2, gi, and g2 are dimensionless. These structure 
functions are related to others in common use by: 



= Pi, W2 



M2 

— P2, Gi 

V 



Go 



M^g2 



V 



;i.i6) 



1.1.4. Scaling and Kinematic Domains 

Our choice of invariant structure functions makes the determination of scaling behavior at 
large almost trivial. In the Bjorken limit where — > cxd and v —>■ 00, x = Q^/2z/ fixed, 
QCD becomes scale invariant up to logarithms of generated by radiative corrections. 
Under a scale transformation, P —>■ AP, q —>■ Xq, and M —>■ AM 7^ M, so a theory with a 
discrete spectrum of massive particles cannot be scale invariant except in a limit in which 
all masses are negligible. Thus no masses can appear in W^" in the Bjorken limit; it must 



7 



be a dimensionless function of P^, q^, S^, and the invariants Q and z/. In particular, it 
cannot depend explicitly on the target mass, M. If, for example, a term like ^"^p" W2 ap- 
peared in W^^'^\ it would violate scale invariance unless W2 vanished like — at large Q^. 
Clearly, the way to avoid such pathological choices of structure functions is to write the 
dimensionless tensor W'^'^ in terms of dimensionless invariant functions using u (or Q^) to 
supply dimensional factors as needed. The immediate conclusion is that the functions Fi, 
F2, gi, and g2 defined in eqs. ( |1.14|) and ( |1.15|) , become functions only of the dimensionless 



ratio X = Q'/2z/, modulo logarithms, in the Bjorken limit, 

Fl(g^^/) ^Fi(x, lnQ2), F2{Q\u) ^ ^2(0;, InQ^) 
^l(Q^^^) ^giix,\nQ^), g2iQ^i^) ^g2{x,\nQ^) 

(1.17) 

as and z/ become large at fixed x. In practice it is observed that for Q'^ > IGeV"^, the struc- 
ture functions depend only very weakly on Q^. Furthermore one observes an approximate 
relationship between Fi and F2, known as the Callan-Gross relation,0 

which indicates that the particles that carry electric charge (the quarks) have spin i. The 
different kinematic domains of interest in inelastic electron scattering are shown in fig. (|^). 

1.1.5. Flavor Generalizations 

Only up, down and strange quarks appear to be important constituents of light hadrons. The 
processes of interest to us, therefore are mediated by currents lying in the SU{3)ji x SU{3)l 
space of u,d,s vector and axial currents, 

(1.19) 

where A*^ for a = 1, 2 ... 8 are the flavor SU{3) matrices, which are normalized to TrA'^A'' = 
26°''^. Note, in particular, that A3 = diag(l, —1, 0) and Ag = diag^(l, 1, —2). In addition one 

has the flavor singlet current J° = \J 2 / 3il)'~^ f/ip , acting like Y^diag(l, 1, 1) in flavor space. 

1.1.6. Cross Section for Electron-Hadron Scattering 

The differential cross section for unpolarized electron-hadron scattering can now be expanded 
in the Lorentz scalar structure functions by contracting the symmetric tensor, eq. ( |1.14|) , with 



the leptonic tensor, eq. ( |1.6D . Likewise the cross section for polarized scattering is obtained 
by contracting the antisymmetric tensor, eq. ( |1.15| ), with the same lepton tensor. The result 
is often quoted in terms of the experimenter's variables, Q', 6, E and E', e.g. for the spin 
average case, 

E'^ { 2Wi{q^, v) sin' - + 1^2(9 , v) cos' - . (1.20) 



dE'dn MQ^ 1 xv-^ ' y 2 ' ' ^ 2 
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elastic scattering Q^=2d 
rescnances2u -Q2=m2 _]y[2 

1 

Bjorken limit 




X 







> f 9 

^ Reggelimit I ^ J^^ 



V 



photoprcxluction 



Fig. 4. Kinematic domains in electron-nucleon scattering. 
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The relative importance of the two terms is difficult to judge. Superficially it looks as though 
Wi and W2 are equally important. On second thought, Wi is multiplied by sin^ | which gets 
small in the Bjorken limit. On third thought, W2 vanishes like To disentangle all this, 
we rewrite da in terms of Fi, F2, x, y, and Q^, where scaling behavior should be manifest. 



dxdyd(P 47r2Q2 U ' 2xy 



(1.21) 



with K = 1 



. No scaling approximations have been made in eq. ( |L2l . Under typical 



experimental conditions y and x are of order unity, though experiments are now being carried 
out at very low-x. Since Fi ~ ^ and F2 ~ const, for small x, the two terms are comparable. 
There is no significant dependence on the azimuthal angle </>, which cannot even be uniquely 
defined for inclusive scattering with an unpolarized target. 

It is clear from the tensor structure of i^u and W^u that no target spin dependent 
effects survive if the beam is unpolarized. Therefore we lose no generality by defining the spin 
dependent cross sectiom Ao" as half the difference between right- and left-handed incident 
electron cross sections,! 



dAa{a) 
dx dy dcf) 



47r2Q2 



cos a 



1 y y ( ^\ 

1 [K — 1 

2 4^ ' 



y 

gi{x, Q^) - -(k - l)g2{x, Q^] 



— sm a cos ( 



[k-1) 1-y 



(k - 1 




gi{x,Q'^) + 5(2(2;, 



;i.22) 



Now the azimuthal angle and the angle, a, between the target spin S and the incident 
electron momentum, k, make non-trivial appearances. These and other kinematic variables 
are defined in fig. (^). Note the following: 

• a is the angle between the spin vector of the target (S) and the incident electron beam 
(/c), not the virtual photon direction (g). 

• is the azimuthal angle between the plane defined by k and k' and the plane defined 
by k and S. 



Eqs. ( |1.21|) and ( p..22| ) are exact (except that lepton masses have been ignored): no 
scaling limit has been taken, k — 1 = M'^Q'^/v'^ = AM'^x'^/Q'^ is a measure of the 
approach to the scaling limit, Q"^ ^ 00. 

To eliminate spin-independent effects one may either (i) subtract cross sections for 
different values of a; (ii) subtract cross sections for right- and left-handed leptons; or 
(iii) measure 0-dependence. 



Notice that effects associated with g2{x, Q'^) are suppressed by a factor V/t — 1 



2Mx 



with 



respect to the dominant structure function gi{x,Q'^). In technical terms, this means that 
effects associated with g2 are "higher twist" — suppressed by a power of Q relative to the 
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Fig. 5. Kinematic variables in polarized lepton scattering from a polarized target. 



Im 




< 




Fig. 6. The optical theorem relates the total cross-section for e^e annihilation to the imaginary part 
of the photon propagator. 

leading phenomena in the Bjorken hmit. However, at 90° the coefficient of the dominant 
term vanishes identically and allows the combination ^gi + g2 to be extracted cleanly at 
large Q^. This is a unique feature of the spin-dependent scattering. Only very rarely, to my 
knowledge, can a higher twist effect be selected by an adroit kinematic arrangement, thereby 
avoiding the difficult process of fitting and subtracting away a leading twist effect to expose 
the higher twist correction underneath. 

1.2. Other Basic Deep Inelastic Processes 
1.2.1. Inclusive e+ e" Annihilation 

In this process an electron with momentum k and a positron with momentum k' annihilate 
to form a massive time- like photon with momentum q = k + k' {Q^ = > 0), which 
decays into an unobserved final state. Through the optical theorem, the total cross section 
is proportional to the imaginary part of the photon propagator (see fig. (j^)). 



where Il{Q^) is the Lorentz scalar spectral function appearing in the photon propagator: 



Usually the data are expressed as a ratio to the pointlike annihilation cross section to muons 






and 
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(to lowest order in aEAd)'- 



Rm ^ ^ ""^^ = i2.iim. (i.26) 

Since the hadronic process is initiated by the creation of a qq pair, R directly measures the 
number of colors. At large it is modified only by perturbative QCD corrections: 

+ quark mass corrections.) (1-27) 



The coefficients in eq. ( |1.27|) are renormalization scheme dependent beyond lowest order. 



Those quoted in eq. ( |1.27| ) were calculated in MS scheme with five flavors." The formula 
for R does not depend on any details of hadronic structure, so it provides an important 
test of QCD (and measurement of agcz))- Similar remarks apply to processes in which 
jets are observed in the final state of e^e~ annihilation. Two jet events have the angular 
distribution that one expects for two spin 1/2 quarks; a third jet is associated with gluonic 
bremsstrahlung. These processes however are not sensitive to the structure of hadrons and 
we will not discuss them further here. 



1.2.2. Inclusive e+e Annihilation with One Observed Hadron 

This process looks very much like a timelike version of deep inelastic scattering. Indeed it 
shares many important characteristics, but it also differs in essential ways. From the point of 
view of a theorist interested in hadron structure, the opportunity to study unstable hadrons 
makes this process very attractive. Deep inelastic scattering from A-hyperons or tt or p- 
mesons will never be more than a gedanken experiment. However, these and other unstable 
hadrons have already been studied in e"'"e~-annihilation. The physical basis of "fragmenta- 
tion" — the process by which a quark created by the current from the vacuum fragments into 
the observed hadron — is not as well understood as DIS, making this an area of considerable 
interest at the present time. 

The kinematics for e~^e~ — > P + X is illustrated in fig. (|^. Once again two kinematic 
invariants, and u = P ■ q, define the process. The limit of interest is Q^, v —>■ oo, with 
z = ^ fixed. The momentum of the virtual photon is time-like, and that makes a major 




Fig. 7. Kinematics for single particle inclusive annihilation — e^e PX. 
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difference as we shall see in §2. The invariants are often expressed in terms of quantities 
measured in the e+e~ center of mass: 

P-q = v = E^ 

2P-q E 

0<z ^ ^ = — <1, 
Ee 

(1.28) 

where E is the energy of the observed hadron. We shall usually be interested the polarization 
dependence, but here we illustrate the kinematics for the simpler, spin-averaged case. The 
cross section can be written as the product of a leptonic l^i, and a hadronic tensor W'^'^: 

The hadronic tensor is determined by the electromagnetic current and depends on two in- 
variant "fragmentation functions" due to current conservation and C, P and T invariance: 

^'"^ = ^T.i'^7iy6\P + Px-q){0\J,\PX),^,^^,{PX\J,\0) 



i- J d'^e^^<Y.iW0\PX)outout{PX\M0m 

F,iz, q-) + 1 (P, - ^ J U - ^gJ F,(., q% 



;i.30) 



In contrast to DIS, the sum over unobserved hadrons X cannot be completed because the 
state |PX)out depends non-trivially on the observed hadron. Even if P and X did not interact, 
Bose or Fermi statistics prevents the states X from being complete. In practice P and X 
interact dynamically, as indicated by the subscript "out". For simplicity we will generally 
suppress this subscript. Thus, e+e~ P+X is not controlled by the product of two operators 
(electroweak currents), a feature which complicates the study of e^e~ P + X significantly. 

If g2 — > cxo at fixed z, the structure functions, Fi and F2 scale (up to logarithmic 
corrections) and obey a "Callan-Gross" relation, Fi + {z/2)F2 ~ l/lng^. In this limit the 
cross section is: 

^ = + (1.31) 

In leading logarithmic order, using "Callan-Gross", the inclusive spectrum reduces to 

1 dcr 2 " , r,, , , 

-— ~-2;Fi(z,g2) + ..., (1.32) 

a dz K 



where R is defined by eq. (|1.26| ). 
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Fig. 8. Kinematics of the Drell-Yan process. 



1.2.3. Lepton Pair Production 

The final process we will consider in detail is massive lepton-pair creation in hadron-hadron 
collisions - the so-called "Drell-Yan" process. The opportunities for study of novel aspects of 
hadron structure by means of polarized Drell-Yan experiments have motivated a major spin 
physics program at RHIC.0 The kinematics of the lepton pair production are illustrated in 
fig. d^). Two hadrons with momenta P and P' collide at a center of mass energy s = {P + 
P'Y = ^E'^m- Two leptons with momenta and fcg respectively are produced. They result 
from the decay of a timelike photon, W^, or carrying a momentum g, with = > 0. 
Two dimensionless scaling variables are defined by x = and y = ^^r^- It is easy to see 
that xy = Q'^/s. The differential cross section is 

da oc LW-"-^-^-^?^. (1.33) 

The decay of the virtual gauge boson is described by the leptonic tensor L^^^, whereas all 
information about the hadronic process are contained in W'^'^: 



W^^ = lsJ2{27ry6\P + P'-q-X) ,^{PP'\J^mX){X\rmPP')-,n 
^ X 

(1.34) 

where the in-state label on \PP') will usually be suppressed. W'^'^ contains many Lorentz 
invariant structure functions Wk- Depending on the experimental circumstances different 
combinations of the Wk and differential cross sections are of interest. As an example we 
consider the inclusive cross section where the lepton momenta have been integrated out, 
leaving da/dq"^, 

^ - gij^ (-t^';) . (1-35) 
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The scaling limit (s, —>■ oo but r = Q^/ s fixed) once again yields a function of the 
dimensionless variables {x and y) modulo logarithms induced by QCD radiative corrections, 
and in this case W,^ is of interest. 

Wi:^W{x,y,\nQ'). (1.36) 

2. Deep Inelastic Processes from a Coordinate Space Viewpoint 

Traditional introductions to the parton model stay fixed in momentum space, where 
they use the device of the "infinite momentum frame" to simplify dynamical arguments. 
More sophistication is necessary to handle the complexities introduced by spin dependence 
and the subdominant effects associated with transverse spin in DIS. It is particularly useful 
to employ coordinate space methods, mixing partoiiDhenomenology with somewhat more 
formal methods of the operator product expansion.li^ Certainly, sophisticated momentum 
space methodslli can achieve the same results. However, it is particularly easy to distin- 
guish and catalogue dominant and sub-dominant contributions using the operator product 
expansion in coordinate space. 

In this section we will explore the coordinate space structure of the hardprocesses 
introduced in §1. Much of this material is to be found in modern field theory texts,Bj however 
there is an advantage to providing a brief, self-contained introduction which stresses only 
those elementary aspects of the formalism that are useful in characterizing deep inelastic 
spin physics. 

2.1. e~^e~ hadrons - The Short- distance Expansion 

Inclusive e^e~ annihilation into hadrons is the simplest process to analyze and illus- 
trates the importance of Wilson's short distance expansion. As shown in §1, this process is 
described by the vacuum expectation value of a current commutator, 

n(Q^)oci^/rf^ee^'^-Mo|[J,(0,^'^(o)]|o). (2.1) 

In the center of mass system we have q = {^/Q^, ). Since the commutator is causal, 

[UO,JM] = for e<0, (2.2) 
then 1^1 < in the integral. Using the symmetry of the commutator one obtains: 

Jo J\^\<^o 

In the high energy limit, Q — oo, sinQ^^ oscillates rapidly, averaging out contributions 
except at the = boundary of the integration region. This argument can be made more 
formal, leading to the conclusion that ~ ^ gives the dominant contribution to the integral. 

Since > |^| we can conclude that e^e~ annihilation into hadrons at high is dominated 
by interactions at short distances, —>■ 0. This is, of course, a Lorentz invariant condition. 

The leading contribution to the annihilation process can now be found via the operator 
product expansion (OPE).^ First postulated by Wilson, the existence of the OPE has been 
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demonstrated to all orders in perturbation theory in renormalizable theories and also in 
various toy models which can be solved exactly. According to the OPE, a product of local 
operators A{^) and -B(O) at short distances (here — >• 0) can be expanded in a series of 
non-singular local operators multiplying c-number singular functions, 

i(O5(0)~5:q,](O^H(0) as e^O- (2.4) 

[a] 

In general the product AB is singular as ^ ^ 0. The substance of the expansion is that 
the singularities can be isolated in the c-number "Wilson coefficients", C[a]- The operators 
in eq. ( p.4| ) are cutoff independent renormalized operators and the Wilson coefficients are 
likewise cutoff independent. 

The behavior of the Wilson coefficients at — follows from dimensional analysis. In 
natural units, all quantities are measured in dimensions of mass to the appropriate power. 
For simplicity, if a quantity, 6 has units m'^*, we write dg = [6]. This is a simple concept, 
not to be confused with more subtle ones like anomalous dimensions or scale dimensions." 
The dimension of all operators of interest to us can be deduced from the fact that charge 
and action are dimensionless. Thus [J'^] = 3 because / d^xJ^{x) = Q. For the quark field 
[ip] = I because the free Dirac action is / d'^xtpi'j ■ dip + . . ., likewise for the gluon field 
strength = 2. Since we normalize our states covariantly, {P\P') = 2E{27i'Y6^{P — P'), 
[{P\P')] = —2. For the vectors [P'^] = [S^] = 1. We see that W^^, is dimensionless, as reflected 
in the form of eqs. ( |1.14| ) and ( |1.15| ). 



Dimensional consistency applied to the OPE requires, 

[i] + [B] = [q„]] + [y. (2.5) 

What can account for the dimensions of the singular function If the operators A and 
13 are finite in the mquark limit, then powers of mquark can only appear in the numerator 
of C. The renormalization scale, /i, necessary to render the theory finite can only appear 
in logarithms (of the form ln(/i^)) order by order in perturbation theory. This leaves the 
coordinate ^ itself to absorb the dimensions. 

^N(0- ^,^Jj_,,^^/ ln--(/.0 + ...) as (2.6) 

The exponent 'jg is the "anomalous dimension" of the operator 6 generated by radiative 
corrections. Without minimizing the importance of these logarithms, we will usually ignore 
them and focus on the gross, power law, behavior required by dimensional analysis. For given 
operators A and i? the leading contribution at short distances comes from that term in the 
OPE having the lowest operator dimension [6[a]]- 

This can now easily be applied to e^e~ annihilation. The dimension of the hadronic 
electromagnetic current is [J^] = 3. No fields have negative dimensions, so the lowest dimen- 
sion operator is the unit operator, 9o = 1, with [1] = 0. The Co(0 ~ 1/^^ aiid the dominant 
contribution in the OPE is, 

(0| [J^(0, J^m |0) ~ ^, modulo logarithms. (2.7) 
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Consequently the current correlation function scales like 

m')-^Jd'^e^^<^^^l, (2.8) 
again modulo logarithms, and the cross section eq. (|1.23|) scales like: 

a (e~^e~ —>■ hadronsj ~ — ^. (2.9) 



The logarithms can be gathered together into powers of ^ as anticipated in eq. ( |1.2 
Of course, having made no attempt to derive the OPE or to study the effects of radiative 
corrections and renormalization in detail, the example of the total e^e~ annihilation cross 
section becomes rather trivial. Nevertheless it provides a useful introduction to the more 
complicated cases which follow. 

2.2. Ip IX - The Light-Cone Expansion 

Next we turn to deep inelastic scattering, which is characterized by two large invariants 
- and u. As we shall see, such processes are dominated by physics close to the hght-cone. 

2.2.1. Light-Cone Coordinates and Formulation of Deep Inelastic Scattering 

The four-momenta and can be used to define a frame and a spatial direction. Without 
loss of generality we can choose our frame such that and q^^ have components only in the 
time and 63 directions. It is helpful to introduce the light-like vectors 

p^ = ^(1,0,0,1), 

:i, 0,0,-1) 



^/2p 

(2.10) 

with = p2 = and n ■ p = 1. Up to the scale factor p, the vectors p'^ and n'^ function as 
unit vectors along opposite tangents to the light-cone. They may be used to expand and 





= ± 




qf" 




l^z/ - ^ 




= + 


M2 


n 

2 



Z/2 + M2Q2 j + ^^2 + JVf2Q2 J (2.11) 



In the Bjorken limit simplifies to 



(2.12) 



limg'^ ~ ( z/ + -M^x ) n^" - xp^ + O { — ] . (2.13) 



1 



p selects a specific frame. For example p = M/ \/2 yields the target rest frame, while 
p ^ oo selects the infinite momentum frame. The decomposition along p^^ and is equivalent 
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to the use of light-cone coordinates, which are defined as follows. An arbitrary four-vector 
a'^ = (a°, a},a^,a?) can be rewritten in terms of the four components = -^{oP ± a^), and 

a-*" = (a^,a^). In this basis, the metric gfj_u has non-zero components, = g |- = 1 and 

9ij ~ so a ■ b = a^hr + a~b^ — ■ b-^. The transformation to light-cone components 

can be recast as an expansion in the basis vectors and n'^, 



a 



V 



+ (v^a+p) + a^^. (2.14) 



With these preliminaries it is easy to find the space-time region which dominates the 
DIS. Consider the hadronic tensor W^"^ defined in eq. (|1.8|) : 

W^"" = ^ f d^^e"'-^ (P| [J'^(0, J"(0)] \P) , (2.15) 
47r J 

Take the Bjorken limit by keeping P fixed and g — > oo. Define 

= r^p^ + An^ + e^^, (2.16) 

we find in the Bjorken limit: 

limq ■ ^ = rjv — x\. (2-17) 

Arguments similar to those used in the previous section show that the integral in eq. ( p.l5| ) 
is dominated bY?7~l/z/~0 and A ~ 1/x, which is equivalent to ^~ ~ and E,^ ~ 1/xp 
respectively.iiS As in the previous case the commutator in eq. ( p.l5|) vanishes unless = 
2\ri — > because of causality. Combining these results we find that the Bjorken limit 
of DIS probes a current correlation function near the light-cone = 0, extending out to 
distances (^^ and of order 

2.2.2. Deep Inelastic Scattering and the Short Distance Expansion 

QCD simplifies at short distances on account of asymptotic freedom. The analysis of e^e~ — *• 
hadrons simplifies greatly for this reason. Deep inelastic scattering is not a short distance 
process; it is light-cone dominated. Nevertheless it can be related to the OPE and to short 
distances with considerable resulting simplification. 

To show this we consider the so-called the Bjorken- Johnson-Low limit (lim5ji).@ii 
This is a somewhat old fashioned method, mostly supplanted by Wilson's operator product 
expansion. It has the virtue that the connection between measurable structure functions and 
local operators is extremely clear (via dispersion relations). Use of the BJL limit prevents 
one making mistakes in subtle cases.ii'0 In the BJL limit one takes g = and g° ioo, 
which yields —>■ —oo and x —>■ —ioo. In the physical region x is restricted to be real and 
between and 1. So the hadronic tensor W^u cannot be measured in the BJL limit. It is 
useful because 1) it is dominated by short distances, and 2) it can be related to Wfj,^, in the 
physical region through dispersion relations. Remember that Wf^u is the imaginary part of 
the forward, virtual Compton amplitude, T^^, by the optical theorem, 

T,.{q\u)=t J d'^e'^<{P\T{J,{OJum\P), (2.18) 
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Im (0 



0) 




For simplicity we suppress Lorentz indices and spin degrees of freedom for a while. Standard 
dispersion theory arguments show that T(g^, u) is an analytic function of u at fixed spacelike 
with branch points on the real-i^ axis at u = the threshold for the elastic process 

7*p — > 7*p. In fig. one can see the physical region of this process and the area of the BJL 
limit in the complex u = 1/x plane. The physical cuts lie on the real axis from - = ±1 to 
±cxo. This means that T(g^, -) is analytic within the unit circle about the origin. The BJL 
limit takes - to zero along the imaginary axis. Thus T^q"^, u) can be expanded in a Taylor 
series in (^) about the origin in the BJL-limit. The coefficients in the Taylor expansion can be 
obtained from the dispersion relation obeyed by T. First remember that the optical theorem 
relates the imaginary part of T to the hadronic tensor W{q'^, v) in the physical region, 

ImT(g^ v) = 47iW{q\ v). (2.19) 

Dispersion theory tells us that an analytic function can be represented in terms of its singu- 
larities in the complex plane,0 in this case the physical cut on the real axis, 

nfM^iTi^'j'^^^. (2.20) 
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Crossing, i.e. T(uj) = T{—uj), has been used. Since T(g^, x) is analytic for |i| < 1 it may be 
expanded in a Taylor series in powers of -: 



BJL 

with 



hmT{q\x)=A ^ M^q') (2.21) 



M"(g2)= f' dxx''-'W{q^,x). (2.22) 







Now consider where the BJL limit leads us in coordinate space. With q = and 
q^ ioo, the factor e*''^ in eq. (|2.18| ) reduces to e"''' and forces to zero. Although 



the time ordered product does not vanish outside the light-cone, it can be exchanged for a 
"retarded commutator" ,1111 which does. Thus ^ forces —>■ and we conclude that 
the BJL-limit takes us to short distances where Wilson's operator product expansion may be 
used. The OPE analysis of the product of currents yields a power series in - multiplying the 
matrix elements of local operators. Identifying terms in this Taylor series with the terms in 



eq. p.21j we obtain the celebrated "moment sum rules" relating integrals over deep inelastic 
structure functions to target matrix elements of local operators. We will not pursue this 
direction further here — it is treated in standard references.00 

2.3. e^e^ — > hX- Once Again, the Light-Cone 

Like deep inelastic scattering, single particle inclusive production in e~^e~ annihilation 
is dominated by the light-cone. However, the operator product expansion does not apply and 
no short distance analysis exists. The process is described by the tensor introduced in §1, 

W,. = ^J d'^e^'^<Y.(^\UO\PX) {PX\jmO) (2.23) 

Once again, the nucleon and photon momenta may be expanded in terms of the light-like 
vectors introduced in §|2.2|, 



M2 

^ 2 
1 



and in the Bjorken limit (Q^, u ^ oo with z finite). 



(2.24) 



/ 1M2\ 1 
limg'' = I z/ - -— j + -p'^. (2.25) 

It is traditional to use the photon rest frame (g = 0,p ~ to analyze the process. However, 
the label on the state, {P^) changes as the limit ^ oo is taken in this frame, making it 
difficult to sort out the important regions of the ^-integration. Things are simpler in a frame 
where P is fixed, e.g. the rest frame of the produced hadron, where p = In such a frame. 
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the analysis of the fourier integral in eq. ( |2.23D proceeds exactly in the same way as for the 
electroproduction process of §p.2|. With 



= rip^" + Xn" + (2.26) 

we find in the Bjorken limit 

\imq-^ = r]u--. (2.27) 

So, z/ — > oo implies — > and A ~ 2;, since z is finite. So light-like separations ~ 
dominate again unless unusual variations occur in the matrix elements 

$:(0|J(O|PX)(PX|J(0)|0), (2.28) 

X 

which will not happen in the frame where P is independent of and u. Also the frequencies 
associated with the states in the sum J2 know nothing about and u, and will 

X 

not spoil the argument. For a contrasting situation see the discussion of Drell-Yan in the 
following sub-section. 

One can thus conclude that light-cone distances dominate fragmentation. However, in 
contrast to DIS the OPE cannot be applied here since the observed hadron state, \P) in- 
terferes with the attempt to complete the sum on X. Nevertheless nearly all of the QCD 
phenomenology developed for DIS can be carried over to this case, primarily using momen- 
tum space methods we will not discuss here.0 In §6 we will see that the limitations on 
the coordinate space analysis do not prevent us from analyzing spin, twist and chirality in 
fragmentation. 

2.4. PP - The Drell-Yan Process 

Finally we consider the Drell-Yan process. Here the relevant hadronic tensor is (see 

§1): 

W^"" = ^sj d^^e'"-'^ {PP'\J^{i)r{Q)\PP') . (2.29) 

It is simplest to consider the case where only the dilepton invariant mass distribution da/dq^ 
is measured, though other observables behave similarly. depends only on Wj^. We define 
a function, W{s, Q^) by integrating PF^ over all with = and g° > 0, 

W{s,Q') = -^J^dU5{q'-Q'){-g^''W,Mq'), (2.30) 

= -ATrh I d^q 6{q^ - Q^) Y,{2tx*)6\P + P' - q - X) 

X {PP'\J^'\X) {X\Jf,\PP') . (2.31) 

The virtual photon's momentum is integrated over all values consistent with the constraint 
q^ = and conservation of energy, ^/Q^ < q^ < {s + (5^)/2v^, which defines the region R. 
If we introduce the function 

Ki^, Q') = e-«-«5(g2 - qyiqo), (2.32) 
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then W{s, Q^) can be written as 

wis,Q') = -s I d'^A^i^,Q') {pnno-m\pn ■ (2.33) 

In the scahng hmit (Q^, s ^ oo, r = Q'^/s fixed) A:^(^, Q^) approaches a well-studied func- 
tion of quantum field theory, the free field singular function, A_|_(^, 



+ ^^^^^(-nK.[Q^e). (2.34) 

A+(^,Q^) is singular on the light-cone and would select out light-cone contributions were 
it not for high frequency variations in the matrix elements. These can occur because the 
hadron momenta P and P' cannot be kept fixed as ^ oo (s ~ 2P ■ P' > Q"^) in any 
frame. Even in free field theory or the parton model the matrix element behaves like 

{PP'\J^{i)JM\PP') ~ / dadl3e'''''<+'^'"<f{a,P), (2.35) 

where /(a, j3) labels the momentum components of the partons that contribute to the current. 
To see that such variation can lead to contributions off the light-cone, consider a frame defined 
through the two vectors 

p'^ = p (1,0, 0,1), 



The hadron and photon momenta can be written as 

s 

S 

gfj. ^ yp^' + xp'^, for Ig-rl < vQ^- 



(2.36) 



(2.37) 



With = rjpf^ + Xp'f^ + the hadronic tensor eq. ( |2.29| ) is then equal to 



Wl^ = J dXdrjd^^T J dad(5 f{a,f5)e'-2^^''-y^^+^^-''^'^\ (2.38) 

Therefore the phases will cancel and the Drell-Yan process will escape from the light-cone if 
a y = 2q ■ P' / s and /3 ~ a; = 2g ■ P/ s. In §5 we will return to this process and see that 
such phases are generated in a natural way. 
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2.5. Dominant and Suhdominant Diagrams 



Guided by our understanding of the regions of coordinate space important for various 
deep inelastic processes, we can return to the more famihar world of Feynman graphs and 
learn which diagrams are likely to give dominant and subdominant contributions. The quarks 
that couple to electroweak currents propagate according to Sf{C), the Feynman propagator. 
In coordinate space, SfH) behaves like at short distances (note Sf{^) ~ / d'^P ^^Jl^ )-, 

• > • ~ Sf (I) ~ ^ . 

C ^ 

Interactions will not increase the singularity. For example, coupling a gluon to the propagat- 
ing quark gives. 




Generally speaking in renormalizable field theories, interactions on propagating lines do not 
increase the order of the short- distance or light-cone singularity by more than logarithmic 
terms beyond free field theory. This can be used as a guideline to estimate the importance 
of different perturbative diagrams for hard processes. 

As a first example, consider e~^e~ hadrons. The total cross section is proportional to 
the vacuum polarization of the photon propagator, whose leading contribution results from 

the quark-antiquark loop fourier transformed, this 0{-^) behavior generates a cross section 
which scales like Radiative corrections introduce logarithmic dependence on In^^/i^, 
where /i^ is the renormalization point, but they do not change the power of the singular- 
ity in a renormalizable theory. The renormalization group may be used to sum classes of 
diagrams giving modifications of the ^ behavior which go like powers of logarithms in an 
asymptotically free theory like QCD. 

In deep inelastic scattering the leading contribution to the cross section or the forward 
Compton amplitude is shown in fig. (0). It dominates because the free quark propagator has 
the greatest possible light-cone singularity. The modifications shown in figs. ([lTl - |l4|) introduce 
only logarithmic modifications of the ~ 1/^'^ singularity. Renormalization group summation 
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Fig. 10. Leading diagram in deep inelastic scattering. The quark propagator between the two currents 
carries the large momentum q^^ and leads to a 1/^^ behavior at small distances. 





Fig. 11. Radiative corrections to the quark propagator lead to O (as(^^)) corrections to the coefficient 
of the leading 1/^'^ term. 



of leading In^^ dependence leads to powers of In^^/i^ but no change in the fundamental 
power singularity. All radiative corrections can be classified in the fashion outlined by figs. 
(01" li)- -^^^ single particle inclusive e~^e~ annihilation one finds in analogy to deep inelastic 
scattering the leading diagram fig. (^) which has a singularity 1/^^ from a propagating 
quark. Radiative corrections can be treated as before. If light-cone dominance were the only 
consideration, the diagram in fig. ([T6|) would dominate the Drell-Yan process. However, if 
one studies the flow of hard momentum, this diagram turns out to be suppressed. The quark 
which brehmsstrahlungs the massive photon must be far off-shell, which is unnatural in 
a hadron-hadron collision. In coordinate space this is reflected by the fact that no large 
phases are generated by the matrix element. The dominant contribution to the Drell-Yan 
process is shown in fig. (|l^). The enclosed parts appear to be identical to the quark-hadron 
amplitude that occurs in the diagram that dominates in deep inelastic scattering. This means 
that at tree level, the same structure functions that appear in deep inelastic scattering also 
contribute to the Drell-Yan process. The diagram should still be dressed with QCD radiative 
corrections. The factorization theorem of QCdEI says that this correspondence survives even 
in the presence of radiative corrections. 
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k2~Q2 
k2>^2 





Toinnnn 



(b) 




Fig. 12. The higher order contributions in the upper part of the diagram (a), where the quark virtuality 
is greater than /x^ lead to more O (q;s(^^) In/x^^^) corrections. The lower radiative corrections (5) can 
be absorbed into the quark-hadron amplitude which will then depend on the renormalization scale fj?. 



k2~Q2 
k2>^2 




k2<^2 




Fig. 13. The upper part (a) here generates a c-number coefficient function for gluonic operators in the 
product of two currents, beginning at O (as(^^) In /x^^^), while the lower part (6) can be absorbed into 
a new gluon-hadron amplitude. 
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A <=>{ <=>{ 

2 J Z 





Fig. 14. These corrections are either gauge artifacts or modifications of lower order in 1/^. 




Fig. 15. The dominant diagram for the single particle inclusive e"^e annihilation. 




Fig. 16. The diagram that has the leading light-cone singularity of the Drell-Yan process, but does not 
dominate. 
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Fig. 17. Dominant diagram of the Drell-Yan process. 

A subtlety of the Drell-Yan process is that the term most singular on the light-cone 
does not dominate, nevertheless the diagram gets its dominant contribution from ^ 0. 
Returning to the definition of W{Q'^,s), we see that Ar{C,,Q^) forces to zero, but the 
phase factors generated by the two separate quark-hadron amplitudes select tangent planes 
to the light-cone that contribute to the ^+ and ^~ integrals respectively. 

One can now generalize these results to other processes, appealing to factorization.!! 
For example fig. (p]8| ) shows the dominant contribution to one particle inclusive deep inelastic 
scattering in the current fragmentation region. Building blocks of these calculations are (a) 
and (6), defined and measured in e~^e~ — > hadrons and DIS respectively. Factorization allows 
them to be carried from one process to another. They are the fundamental objects of study 
in hard inclusive QCD and command our attention. 

3. Deep Inelastic Scattering and Generalized Distribution Functions I 

In this first of two sections on deep inelastic scattering the focus will be on developing 
the tools necessary to perform a complete classification of effects at leading and next-to- 
leading order in We begin with some simple considerations of dimensional analysis, which 
we then apply to introduce the operator product expansion (OPE) and introduce the concept 
of "twist" which is useful to classify contributions to hard processes. To proceed further we 
must understand how to treat the Dirac structure of quark fields on the light-cone. This leads 
us briefly to explore light-cone quantization and introduce helicity, chirality and transversity 
as they apply to this problem. We will then look in some detail at a typical leading twist 
(C(l)) and next-to-leading twist distribution before attacking the complete problem 

in 54. 
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Fig. 18. Dominant diagram for the one particle inclusive deep inelastic scattering. 
3.1. Twist 

In §2 we introduced the operator product expansion (OPE) as a tool for analyzing 
e'^e~ — > hadrons, where the operator with lowest dimension dominates. We also argued that 
light-like distances (^^ ~ 0) dominate deep inelastic scattering. However, operators of high 
dimension can be important in this case. Instead a new quantum number, "twist", related 
to both the dimension and spin of an operator, orders the dominant effects. 



3.1.1. Twist and the OPE 

As we learned in §2, the hadronic structure tensor of deep inelastic scattering, 

AnW,, = / d'^e^^-^P, S\[J,{0: J.(0)]|P, S) (3.1) 

is dominated by .^^ ~ in the — ^ oo limit. To make use of this, we expand the current 
commutator in terms of decreasing singularity around = 0, 

[^(0, ^(0)] ~ E^[^l(e')e^ ■ ■ ■e-^^Mi.-.Mn.lO), (3.2) 

m 

where ^^^^...^^^ (0) are local operators, and K[0]{^'^) are singular c- number functions that can be 
ordered according to their degree of singularity at = 0. Operators of the same singularity 
at {2 ^ will be of the same importance as ^^^^ even though numerator factors of 
render some less singular than others as 0. For simplicity we have suppressed all labels, 
including spin, on the currents J{^) . It often convenient (and sometimes essential) to regroup 
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the terms in eq. ( |3.2| ) so that the operators 6^^ ^^^^^ are traceless {i.e. g^^^^O^^,,,^^^ = 0, etc.) 
and symmetric in their Lorentz indices. We will assume this has been done. 
Substituting the OPE into the definition of the structure function gives: 

AnW = I rfV^-^E^[^l(e')r • • {P\e,,...,.,mP). (3.3) 

m 

where the matrix elements have the form: 

{P\0,....,r.,mP) = P,. ■ ■■P.^^M'^-^o-'fe + .... (3.4) 

The . . . represent several types of terms which are less important in the Bjorken limit. We 
will return to them after looking at the dominant term. 

Note that the power of a mass scale which appears in this expression is determined by 
dimensional analysis alone. We use the parameter M generically for a typical hadronic mass 
scale M ~ Aqcd ~ Relg ~ Mn/3. The power with which M occurs defines the twist of the 
operator 9, 

te = dg- ng (3.5) 

The degree of the light-cone singularity of K^g] ~ ^~^+te jg ^igQ determined by dimensional 
analysis and depends only on the twist, tg. 

To carry out the fourier transformation make the substitution, 



which yields. 



e - -2.g'^^ (3.6) 



^-^-Ei^T'^izTfe (3.7) 



{e} 

,2 



X 



So the importance of an operator as cxd is determined by its twist. As we shall see, it 

is typical for towers of operators with the same twist (and other quantum numbers such as 
flavor) and increasing spin to appear in the OPE. Then it is convenient to sum over spin — 
J2ne fe-^" = fe{.x) - where we now use the label 6 to refer to the entire tower of operators. 

The effect of radiative corrections is to introduce logarithmic dependence on into the 
function fg. Note however that the power law dependence on Q"^ is fixed by twist through 
dimensional analysis. Let us now return to the terms omitted in eq. ( p.4| ). These include 
terms like P^^ . . . -P^„_25'/i„_i/^„M^ that make the expression traceless. It is easy to see that 

these contribute at most corrections of order to the term we have kept. To carry through 
a complete analysis beyond order ^ it is necessary to keep careful track of these terms. This, 
and the fact that interesting spin effects appear at C'(^), are the reasons we do not consider 
O(^) here. 

The lowest twist operator towers in QCD have tg = 2 and scale - modulo logarithms - 
in the Bjorken limit. This reflects the underlying scale invariance of the classical lagrangian. 
The matrix elements of higher twist operators, or the higher twist manifestations of twist-two 
operators are invariably signalled by the appearance of positive powers of mass in expressions 
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analogous to eq. ( ^.41 ) . Dimensional analysis then forces compensating factors of large kine- 
matic invariants in the denominator, suppressing the contribution. The simple conclusion 
is that we can order the importance of effects in the deep inelastic limit simply by keeping 
track of masses we are forced to introduce into the numerators of parton-hadron amplitudes 
in order to maintain the correct dimensions. 



3.1.2. Examples and a Working Redefinition of Twist 

To make the preceding discussion clearer, here are some explicit examples from free field 
theories. These examples are not only pedagogical - the second one generates the leading 
twist effects in QCD up to logarithms. The light-cone singularities can be isolated easily. For 
the time ordered product of two scalar currents built from scalar fields, J{S,) =: 0(^)</)(^) :, 
one can use Wick's theorem to show that, 

T( J(0 J(0)) = -2A^(0 + 4zA^(0 : 0(00(0) : + : 0(00(00(0)0(0) : , (3.8) 

where the normal ordering operation is sufficient to render the operator products finite (in 
free field theory) as ^ ^0, and 

for a massless scalar field. To finally obtain the form of eq. ( p.2|) , simply Taylor expand the 
bilocal operators - 

: 0(00(0) := E • • • U ■■ id'' ■ ■ ■ d^-m}m ■■ (3.10) 

n 

The current associated with a vector flavor symmetry of a fermion field is 

^,"(0=:^(0y7.^(0:- (3.11) 

Making use of the identity 

[V'lV'l, ^20'2] = ^2}^2 - ^2{0'2, ^l}^l (3.12) 

(because {ipi^ ^/'2} = in free field theory), and, for a massless field, 

mim) = ^^<^om% (3.13) 

one can now express the commutator of two currents in terms of bilocal operators:ii 

r\Q)\ = -j- {d,e{^o)S{e)) [d^'^'S^^'-'-Ali^, 0) - td^'^^e^^'^'-Sl^^, 0) 

+ if''^'S>""'''S'^{^, 0) + f''^''e>"""'Al'{^, 0)1 , (3.14) 
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where the Lorentz structure is spht into a symmetric and an antisymmetric part according 
to: 



-Tr7^7p7^7<^ = 9^Lp9ua + QpaQvp - 9^Lu9ap, (3.15) 
and the flavor structure is spht in a similar way: 

A'^A'' = [d^^'X' + if^^X") . (3.16) 
The symmetric and anti-symmetric vector and axial bilocal currents are deflned by, 

s/Af^ ^ mY^ab']m ± mY^ab'm) (3.17) 



Once again the form of eq. (|3.2| ) is obtained by Taylor expanding the bilocal operators. 

We have presented these formulas in their full complexity because they summarize the 
algebra of free quarks at short distances. All of the traditional results of the quark parton 
model applied to DIS (scaling relations, the Adler, Bjorken, Gross-Llewellyn Smith and other 
sum rules, the Callan Gross relation, etc.) can be obtained directly from these relations.^ 

The steps of flrst expanding the bilocal operators, then resumming the tower after 
fourier transformation are very inefficient. Clearly it should be possible to work directly with 
the bilocal operators. The twist content of a bilocal operator is somewhat more complicated 
than that of a local operator. Consider, for example, the bilocal current, ip{0)'f'^ip{C,), which 
occurs in eq. (|3.14|) . The operator has dimension three and, were it a local operator, it would 



have spin-one. In fact it sums an infinite tower of operators of increasing spin and dimension, 
with t >2. For example at short distance one can write: 

mri^io = mrm+^umrd''m + --- (3.18) 

= J'^(0)+e,r"(0) + .... (3.19) 

J'^(O) is traceless, symmetric and local and has twist-two. The operator 6^'^ can be decom- 
posed into a traceless operator and a "trace" : 

9^"" = 1^9^" - ^g'"'9^^^ + ^g^''9l (3.20) 

The first term is traceless, symmetric, with twist-two. The second operator has spin-0, hence 
its twist is four. Further terms in the Taylor expansion of the bilocal operator each yield a 
tower of local operators beginning at twist-two and increasing in steps of two. 

Up to now we have used twist only in the sense in which it was originally introduced — 
tg = dg — ng. In practice, twist is used in a less formal way, to denote the order in ^ (modulo 
logarithms) at which a particular effect is seen in a particular experiment. If it behaves like 
{1/Q'^y, then the object of interest is said to have twist t = 2 + 2p. A traceless symmetric 
operator of twist t will generate contributions that go like (1/Q^)*^^~*\ (l/Q^)*^^"*-* ... as we 
saw explicitly for the operator 9^^- Although the two meanings of twist do not coincide 
perfectly, both are in common use. 
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We will make a definition of the twist of an invariant matrix element of a light-cone 
bilocal operators, that determines the scaling behavior of the matrix element. Matrix ele- 
ments of operators like e.g. ip{0)ji^ip{Xn) are the basic building blocks of the description of 
hard processes in QCD. So we will call "twist" the order in M/Q at which an operator ma- 
trix element contributes to deep inelastic processes. A few virtues of our working definition 
are a) that it is easily read off by inspection of matrix elements; b) that it directly corre- 
sponds to suppression in hard processes; and c) that effects we label twist-t never enter hard 
processes with suppression less than (M/Q)*"^. The twist we associate with the invariant 
matrix element of a specific bilocal operator can be determined simply by considering the 
powers of mass which must be introduced to perform a Lorentz-tensor decomposition of the 
matrix element. The powers of mass carry through the entire calculation to the end where 
each power is compensated by a power of Q in the denominator. Twist-two results in no 
suppression, therefore t — 2 is to be associated with the number of powers of mass introduced 
in the tensor decomposition of a matrix element. 

The method is best explained by example. Consider the spin averaged matrix element 
of the bilocal current, ■0(An)7'*'0(O) on the light-cone, 

{P\ij{Xn)rm\P) =P^h{\)+n^M^h{\), (3.21) 

where the factor of must be introduced because [n^] = — 1. The twist of the first term is 
two but, due to the appearance of the factor M^, the twist of the second term is four. In a 
physical application we assert that the factor of will survive all manipulations and appear 
in the result compensated dimensionally by a factor of l/Q^. Note that it is possible for /i 
to pick up multiplicative factors of jQ^ during a calculation. Twist tells us the leading, 
not the exclusive, dependence of an invariant piece of a light-cone bilocal operator. As a 
second example, consider 

(P|V'(An)^(0)|P) = Me(A). (3.22) 

e(A) has twist- i/iree due to the factor M which must be introduced to preserve dimensions. 
Finally, consider the matrix element of a gluonic operator 

(P|G'/(An)G'«,(0)|P) = p^p./i(A) + (p^n, +p,n^)/2(A)M2 + n^n,h{\)M\ (3.23) 

which has a twist content that can be worked out by the reader. 

3.1.3. Spin and Twist 

Counting twist in the case of polarized targets (or fragments) has an added complication. The 
Lorentz tensors which describe a hadron's spin can appear in the Lorentz decomposition of 
matrix elements — their role in determining twist must be explained. The objects of interest 
in polarized scattering (or fragmentation) are forward scattering matrix elements on a null 
plane: (P, e|^(An, 0)|P, e). The matrix elements are bilinear in e and e*, where e and e* are the 
generalized spinors describing the target (Dirac spinors for spin 1/2, polarization vectors for 
spin 1, etc}). The matrix element is a tensor function of e and e*. For spin- 1/2 the only (non- 
trivial) tensors which can be built from wxu are u'^^u — 2P^, a vector, and u-y^^^u — 25''*, 



32 



an axial vector. We have already analyzed P'^ (it gets decomposed into p'^ and n^). To expose 
the twist content of terms proportional to S'^, express it in terms of and n'^: 



S^' = {S■ n)p'' + {S ■ p)n>' + S^. (3.24) 

Since S ■ p = —M'^S ■ n/2, it is clear that the second term contributes at twist-four. The 
transverse spin term is more subtle. Because we have chosen to normalize S"^ = — M^, [S±] = 
1 and because there are no transverse momenta in the problem, S'^ contains a hidden factor 
of the target mass. At the end of the day this factor will manifest itself in a suppression by 
So we conclude that appearances of S± accompany twist-three distributions. An example 
is provided by: 

{PS\ij{0)ri5H^n)\PS) = S-n p^g,{X) + St_ grW +S-p n''g^{\). (3.25) 

According to dimensional analysis, gi is a twist-two object, gx has twist-three and g^i is a 
twist-four function. When combined with the analysis of the following section, one finds that 
the function we have labeled gi is the scaling limit of the "(/i" defined in §1. Similarly, gT 
turns out to be gi + g2. We discard g^, because we are not concerned with twist-four. 

The same method of analysis can be extended to higher spins. For a spin-1 target, all 
polarization information is contained in the spin-density matrix rj^y = e^e*, which contains 
scalar (e* ■ e), vector (S^ = ^e^i,„/3P''e*"e^), and tensor (17^^ = e*^e^ + e*e^ + g^,u^) polar- 
ization information. Note [5''^] = 1 and [fj^,^] = 2. To determine the twist of the associated 
distributions, rj must be projected along p'^, n^^ and transverse directions^ For even higher 
spins a multipole analysis is more streamlined. @ 

3.2. Dominant Diagram in Coordinate Space 

As a final, and physically important example, we take the dominant diagram identified 
in §2 and use coordinate space methods to compute it. Since the quark that propagates 
between currents suffers no interactions (we are ignoring gluon radiative corrections here), 
we may use free field theory. Working out the commutator of free currents, we get 

w^t = ^Jd'^e^^<{p,s\[r^{0,Jtm\P,s) 

d'^e"^<d^{6{e)e{e)) {S,,..d^'^PS\A^^{^,0)\PS) + ...}, (3.26) 



which corresponds to the handbag diagram of fig. (p!oD. The . . . represent three more terms 



given in eq. ( |3.14| ). This simple free- field picture is modified by: 

• vertex and self energy corrections, which modify the singular function (fig. (p!T|)). They 
give rise to logarithmic corrections, as do the dominant parts of 



ladder graphs (fig. ([T^) ), and 



box graphs, which mix in gluons at 0{as) (fig. (p^). Finally, 

in order to preserve color gauge invariance, one has to remember that the quark prop- 
agates in a gluon background (fig. (|lj 
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On account of the last point, the singular function of free field theory, {iIj{^),iIj{0)} = 
2^^e(^o)^(^^)) must be changed to 

{^^(0,^^(0)} - ^M^o)S{e)V (^expzJ^^dCA.iO^ , (3.27) 

which is the quark propagator in a background gluon field. [The path ordering {V) is nec- 
essary because A^{() is a matrix in color space.] The color field A'^ is that generated by 
remnants of the target nucleon and must be viewed as an operator sandwiched between the 
target hadron states. The bilocal operators in eq. ( p. 26 ) are therefore replaced by, 

where F stands for whatever color /flavor /Dirac matrices appear between ip and ip. The 
5-function in eq. ( p.2(: ) selects the light-cone. If we expand C,'^ about the null plane, = 
Xn'^ + ^'^, it is easy to see that the terms involving are twist-four and higher. One therefore 
has: 

{PSim-P (exp -t rfC^A^(C)) ^(01^-^) = 

{PS\^p{Q)v{ex-p-i drn- A{Tn)^^l){\n)\PS) + (3.29) 

where the . . . represent the parts that vanish on the light-cone and have a twist > 4. In the 
light-cone gauge n ■ A = 0, explicit reference to gluons disappears. However, the inclusion 
of the "Wilson link", P (exp^/o^dCM^(C)), is essential in generating higher twist {t > 4) 
gluon corrections. 

In the unpolarized case, the twist expansion of the bilocal operator matrix element 

gives 

r fix 

/ e'^-{P\^^{Oh^^a{Xn)\P) = 2ha{x)p, + 2M^fUxW (3-30) 
and, carrying out the fourier transform in eq. ( |3.26|) , we find 

Fi{x) = W^lifii^)'fii'^))^ (3-31) 

where a = u,d,s, ... is the flavor index. The interpretation in terms of the parton model will 
be given later in this section. 

A brief summary to this point is: Up to and including twist-three the basic objects of 
analysis in DIS are forward matrix elements of bilocal products of fields on the light-cone 
and in light- cone gauge, 

r d\ 

r{x) = / ^e*^"(P5|^(0)r^(An)|P5). (3.32) 

Remember, that important InQ^ radiative corrections have been ignored in pursuit of the 
twist and spin dependence. 
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3.3. Learning from Light- Cone Quantization 

Since the dominant contribution to DIS comes from the hght-cone, it is natural to 
consider a dynamical formulation in which the light-cone plays a special role. At the birth of 
deep inelastic physics it was recognized that field theories simplify in some important ways 
if they are quantized "on the light-cone" rather than at equal times. Unfortunately some 
features which are simple at equal times become difficult on the light-cone. Certainly, as we 
shall see, there is much insight to be gained by considering deep inelastic processes using 
light-cone quantization. The larger question - whether QCD simplifies in essential ways when 
quantized on the light-cone - will not be pursued here. 

Field theories may be quantized by imposing canonical equal-time commutation (or 



quires that any other space-like hyperplane in Minkowski space would serve as well as = 0. 
A null-plane, such as ^ ■ n = is the limit of a sequence of space-like surfaces, and includes 
points that are causally connected. Although a field theory quantized on at ^ ■ n = could 
differ from one quantized at = 0, they coincide for all examples of which I am aware. Let 
us study what happens if we attempt to quantize fermions on the surface = O.H First we 
must introduce and familiarize ourselves with the unusual kinematics of the light-cone. 

3.3.1. Light-Cone Kinematics 

We have previously introduced light-cone coordinates = :^(^° ± and ^-^ = {^^,^'^), 
and the metric g^^, with = = 1, and gij = The partially off-diagonal structure 
of g makes raising and lowering indices confusing, viz., a'^ = a^, and so forth. So we work 
with upper (contravariant) indices as much as possible. 

Quantizing at (say) = 0, we are committed to C,^ as our evolution variable (just 
as quantization at = fixes as the "time"). ^~ and are therefore kinematic, not 
dynamical variables. The conjugate momenta p'^ and jr^ parameterize the fourier decom- 
position of the independent light-cone fields, just like p in ordinary quantization. p~ is the 
"Hamiltonian" for light-cone dynamics. 

3.3.2. Dirac Algebra on the Light-Cone 

The usual selection of 7° = diag(l, 1, —1 — 1) is prejudiced toward equal time quantization. 
Then a (anti-) particle at rest has only ("lower") "upper" components in its Dirac spinor. 
Much of our analysis is simplified by choosing a representation for the Dirac matrices tailored 
to the light-cone.Ea To represent Dirac matrices compactly, we use the "bispinor" notation: 
let (a^, cr^, 0"^) and (p^,p^,p^) be two copies of the standard (2 x 2) Pauh matrices. A 4 x 4 
Dirac matrix can be represented as ®aK p controls the upper- versus-lower two-component 
space; a controls the inner two-component space. An example will clarify the notation: the 
Dirac-Pauli representation used, for example, by Bjorken and Drell is. 
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iId = = (J J) 

The light-cone representation useful for us is instead, 



(3.33) 



r = l®iB^=r^ (3.34) 



la 
-a^ 






(3.35) 



where ± = 1 or 2. It is easy to check that eq. ( |3.35|) satisfy the usual algebra, {7^, 7^^} = 2g^'-', 
and 7^ = i7°7"^7^7^. 

Operators which project on the upper and lower two component subspaces play a 
central role in light-cone dynamics. Define V± by, 

'P± = i7V = ^(l±«3) 



with the properties: 



7^ = 7f(^°±^')' 



V_V^ = v+V- = 

vl = v± 

= 1 

+ vo y lo 1 



(3.36) 



(3.37) 



The "light-cone projections" of the Dirac field, ip+ = V+ip and ip- = V^ip are known as the 
"good" and "bad" light-cone components of respectively. To save on subscripts we shall 
frequently replace iIj± as follows, 

^+^</) iJ-^x (3.38) 
3.3.3. Independent Degrees of Freedom 

The importance of V± becomes clear when they are used to project the Dirac equation down 
to two two-component equations, 

= -i^-D^x + mx, (3.39) 
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where D± = — igA^. In the hght-cone gauge A'^ = 0. is the evolution ("time") 
parameter, but the first of eq. ( p.39| ) only involves d/d^~ , so it appears that x is not an 
independent dynamical field. Instead the Dirac equation constrains x in terms of and A± 
at fixed 

d 

i7" = -7^ ■ D^<P + m(f) (3.40) 

The longitudinal component of the electric field in electrodynamic is similarly constrained 
{i.e. determined at any time) by Gauss's Law in Coulomb gauge, V ■ E = p. Study of the 
gluon equations of motion indicates that A~ is also a constrained variable. The independent 
fields are therefore and A±. x should be regarded as composite — as specified by eq. ( p.40|) 
-X = ^[^,Aa_]. 

By the way, the reduction of the four- component Dirac field to two propagating degrees 
of freedom is not unique to light-cone quantization. In the usual treatment of the Dirac 
equation one finds only two solutions for each energy and momentum, corresponding to 
the two spin states of a spin-1/2 particle. The two-degrees of freedom corresponding to 
the antiparticle are found in the solution with energy —E and momentum —p. In fact, the 
Dirac equation in momentum space is literally written in the form of a projection operation, 
A_iIj = 0, where A± = -^[f ^m) projects out two of the four components of the Dirac 
spinor. 

Although the complete quantization of QCD requires much more work, the implication 
for the Dirac field is already clear: the good components should be regarded as independent 
propagating degrees of freedom; the had components are dependent fields - actually quark- 
gluon composites. 

The classification of quark spin states depends on the Dirac matrices which a) commute 
with V± and b) commute with one-another. Returning to eq. (|3.35| ) we see that 7^,7^, and 



75 commute with V±. Furthermore, the component of the generator of spin- rotations along 
the es-direction, 

S^-^[V^f]=fn' (3-41) 



also commutes with V±. Note that for a Dirac particle with momentum in the ea-direction, E'^ 
measures the helicity. This set of operators suggests two different maximal sets of commuting 
observables: 

• Diagonalize 75 and — a chirality or helicity basis, or 

• Diagonalize (or equivalently, 7^) - a transversity basis. 
Let us consider these in turn - 



Helicity Basis In the helicity basis, both the good and bad components of ip carry helicity 
labels - the eigenvalues of S^, 



X4 



(3.42) 
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Note that upper and lower components of ijj correspond to good and bad light-cone com- 
ponents respectively. Referring back to the form of 75, eq. (|3.35|) , we see that helicity and 
chirality are identical for the good components of but opposite for the bad components, 



-<P- 

-x+ 

\+X-J 



(3.43) 



This may look strange at first, but it follows immediately from the composite nature of x- 
A quantum of x with positive helicity is actually a composite of a transverse gluon and a 
quantum of 0. Since the gluon carries helicity-one (but no chirality), angular momentum 
conservation requires that the 0-quantum have negative helicity and therefore negative chi- 
rality. Remembering this association will help sort out the chirality and helicity selection 
rules which appear in the following sections. 

Transversity Basis Alternatively, we can diagonalize one of the transverse 7-matrices, 
to be specific, 7^. We define eigenstates of the transverse spin-projection operators, Q± = 
^{l T 757"*^) 5 (which commute with V±), 

Q+<P = 0± (3.44) 
Q-0 = 0T, (3.45) 

and similarly for x- Of course 0t/± are linear combinations of (l)±. Note however that 0t/± 
are not eigenstates of the transverse spin operator T} = ( ^ ^ *n V which is not diagonal 



. 

in the basis of good and bad components of ip. So we have to be careful that we do not 
carelessly confuse transversity , the quantum number associated with Q±, which is simple in 
this picture, with transverse spin, which is not. 

34. The Parton Model 

Following the path we are on, the parton model is merely the light-cone Fock space 
decomposition of the matrix elements which control hard processes. Since we have both the 
matrix elements and the Fock space in hand, it is straight-forward to construct the parton 
model. We will verify that the parton interpretation emerges as expected for twist-two and 
then explore twist-three. The reader should beware that twist-four is considerably more 
complicated. A parton model picture of twist-four does exist, however much work is required 
to make it obvious.ii'0'0 

The Fock space in the two bases can be constructed by defining operators that create 
the appropriate components of (f). In the helicity basis we define K^{k~^, kx) to create a right- 
handed (positive helicity) component of and V'{k'^,ki_) to create a left-handed (negative 
helicity) component of 0, and R\ and Lj, which do the same for the antiparticle field 0c. In 
the transversity basis we define the operators a^{k^,kx) and I3\k^,kx) that create the _l 
and T components of 0, respectively. 
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Fig. 19. The classical parton model. 



3.4.1. Twist-Two 

We begin witli tlie simplest case - the spin average, twist-two deep inelastic scattering which 
is controlled by the bilocal operator defined in eq. ( |3.3CI| ), 

^e''^PS\miM^^)\PS) = 2 {fi{x)p, + h{x)n^} . (3.46) 
We project out the twist-two part, f\, by contracting with n'^. 

/i(^) = U ^e^'^PS\i'mi'i>^n)\PS) (3.47) 
Z J Zn 

1 r dX 



J -e''^PS\<p\0)<P{Xn)\PS). (3.48) 



where we have used the Dirac algebra to express the quark field in terms of its light-cone 
components. Notice that only the dynamically independent "good" light-cone components 
occur. If we make a momentum (/c"*", k±) decomposition of (p and separate helicity states, we 
find, 

/i(x) = - [ (fk^{P\R\xp, k^)R{xp, k^) + L\xp, k±)L{xp, k^)\P), (3.49) 
X J 

This is the parton model as illustrated in fig. ([T^): /i is expressed as a sum of proba- 
bilities to find a (light-cone quantized) quark with p~^ = xp and any transverse momentum, 
summed over helicities and weighted by the phase space factor 1/x. Perhaps the reader is 
more familiar with the "infinite momentum frame" form of the model, where /i is written as 
the sum of probabilities to find an (equal-time quantized) quark with a fraction x of the tar- 
get's (infinite) longitudinal momentum. The two formulations are equivalent since the boost 
to an infinite momentum frame is equivalent to a light-cone formulation. Since eq. ( p. 491) 



is valid in any frame, it can be used in {e.g.)th.e lab frame to provide parton distributions 
which can be associated with quark models.EJ One must, however, be careful to remember 
that the fields in eq. ( |3.49D are good light-cone Dirac components quantized at equal not 
equal 

An identical calculation for a; < captures antiquark operators and leads to the stan- 
dard crossing relation for /i, 

fiix) = -U-x) (3.50) 

where /i is given by eq. ( p.49| ) with R ^ Rc and L ^ L^. The other (spin-dependent) quark 
distributions are explored in the following Section. 
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3.4.2. Twist-Three 



Now let us apply the same analysis to the simplest twist-three distribution function, 

^e'^^P\i>{0)^{Xn)\P) = 2Me{x). (3.51) 

Decomposing in terms of and x? we find 

eix) = ^J ^e^^-(P|0t(o)7\(An) + xHOh'<P{\n)\P) , (3.52) 

which contains the dynamically dependent operator x- If we use the constraint to eliminate 
X we obtain 

e(x) = / ^e^'^Pim^I? ±{^n)<j>{Xn)\P) +h.c. (3.53) 

So e{x) is really a qnaik-gluon correlation function on the light-cone. It has no simple Fock- 
space interpretation in terms of quarks alone, despite the apparently simple form of eq. ( |3.51| ). 

We have happened upon a general (and very useful) result: Every factor of x the 
light-cone decomposition of a light-cone correlation function contributes an additional unit 
of twist to the associated distribution function, 

Twist-2 
(p^X^ Twist-3 
X^X"^ Twist-4. 

(3.54) 

Likewise each unit of twist introduces an additional independent field in the null plane 
correlator: 

Twist-2 (f)^{0)(j){Xn) 
Twist-3 (f)\0)D^{Xn)(f){fm) 
Twist-4 (f)\0)D±{Xn)D±{i2n)(j){iyn). 

(3.55) 

It is as if had twist-one and x h^-d twist-two. 

Twist-three is tractable, using the methods that have been developed in these lectures. 
Twist-four requires a more extensive analysis based on operator product expansion methods 
developed during the 1980's..ii'Ei With these general tools in hand, we turn in the next 
section to the analysis of the specific distributions which appear in deep inelastic scattering 
of leptons. 

4. Deep Inelastic Scattering and Generalized Distribution Functions II 

In this section we use the tools developed in §3 to classify and interpret the quark 
distribution functions which appear in the analysis of DIS. The topics will include the clas- 
sification and parton interpretation of the three leading twist quark distribution functions; 
a discussion of the physics of the less well known transversity distribution, hi, a review of 
transverse spin in hard processes; a short digression on higher spin targets and gluon distri- 
bution functions; and a summary of the physics associated with the twist-three transverse 
spin distribution, g2. 
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(b) 

Fig. 20. Helicity structure a) of the parton-hadron forward scattering amplitude; b) of the n-channel 
discontinuity which contributes to a parton distribution function. 

4.1. Helicity Amplitudes 

Part of the task is simply to enumerate the independent distribution functions at twist- 
two and twist-three. This is simphfied by viewing distribution functions as discontinuities in 
forward parton- (quark or gluon) hadron scattering. Suppressing all momentum indices, each 
quark distribution can be labeled by four helicities: a target of helicity A emits a parton of 
helicity A which then participates in some hard scattering process. The resulting parton with 
helicity A' is reabsorbed by a hadron of helicity A'. The process of interest to us is actually a 
u-channel discontinuity of the forward parton-hadron scattering amplitude ^aa,A'A' as shown 
in fig. (pOl). Note the ordering of indices - although A and A' are the incoming helicities, it is 
convenient to label the amplitude in the sequence: initial hadron, struck quark, final hadron, 
returned quark. 

Since the parton-hadron amplitude results from squaring something like {X\iIj\PS), 
the amplitude must be diagonal in the target spin. However spin eigenstates (in particular, 
transverse spin eigenstates) are linear superpositions of helicity eigenstates, so the {A} do 
not have to be diagonal in the target helicity. Only forward scattering is of interest, so the 



41 



initial and final helicities must be the same, 

A + A' = A' + A. (4.1) 

Also, the parity and time reversal invariance of the strong interactions place constraints on 
the {A}, 

Aax,a'x' — A-A-x-A'-x' (4.2) 
«4aa,a'a' = «4a'a',aa (4.3) 

respectively. 

Clearly the helicity counting outlined above applies equally well to good and bad light- 
cone components of quark or gluon fields. Therefore we can use it together with the methods 
of the previous section to enumerate quark distribution function through twist-three. To work 
through twist-three we will have to consider the case of one good and one bad light-cone 
component. We will identify any bad light-cone fields in helicity amplitudes by an asterix on 
the heficity label. Thus Aqi qi* corresponds to emission of a good fight-cone component and 
absorption of a bad one. 

4-2. Quark Distributions in Targets with Spin-0, 1/2 and 1 
4.2.1. Spin-0 Target 

Only A = is available. Parity equates «^o|,oi ^o-i,o-i- Time reversal equates Aq}_ qI* 
and Aqi* qi. So there is only one distribution function at twist-two and one at twist-three. 

The twist-two function is none other than /i associated with the bifocal operator ijj{0)7^'ijj{Xn) 
and conserves quark chirafity ("chiral even"). The twist-three function is e, associated with 
the scalar bfiocal operator '0(O)'0(An) and flips quark chirality ("chiral odd"). These prop- 
erties are summarized by 



Twist 


A 


A 


A' 


A' 


Chirality 


Two 





1/2 





1/2 


Even 


Tliroo 





1/2* 





1/2 


Odd 



Table 1. Quark distributions in a spin-0 hadron through twist-three. 
4.2.2. Spin-1/2 Target 

In the spin-1/2 case, for each twist assignment there are three independent helicity ampli- 
tudes. The reader may wish to verify that parity and time reversal invariance relate the many 
helicity amplitudes to the six listed in the table below (through twist-three). We leave the 
interpretation of these six distribution functions for the next section where they are discussed 
in detail. 
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Twist 


A 


A 


A' 


A' 


Chirahty 


Two 


1/2 


1/2 


1/2 


1/2 


Even 


Two 


1/2 


-1/2 


1/2 


-1/2 


Even 


Two 


1/2 


1/2 


-1/2 


-1/2 


Odd 


Three 


1/2 


1/2* 


1/2 


1/2 


Odd 


Three 


1/2 


-1/2* 


1/2 


-1/2 


Odd 


Three 


1/2 


1/2* 


-1/2 


-1/2 


Even 



Table 2. Quark distributions in a spin-1/2 hadron through twist-three. 
4.2.3. Spin-1 Target 

A massive spin-one target has three independent helicity states. A new complication ap- 
pears at twist-three: two helicity flip distributions arise which are not related by any of 
the symmetries of QCD. One can easily check that no such complication occurs for spin |. 
There is much interesting physics in these spin-one structure functions, however time will 
not permit us to work through it here. Instead we refer the interested reader to the original 
literature.iiEi 



Twist 


A 


A 


A' 


A' 


Chirality 


Two 


1 


1/2 


1 


1/2 


Even 


Two 


1 


-1/2 


1 


-1/2 


Even 


Two 





1/2 





1/2 


Even 


Two 





-1/2 


1 


1/2 


Odd 


Three 


1 


1/2* 


1 


1/2 


Even 


Three 


1 


-1/2* 


1 


-1/2 


Even 


Three 





1/2* 





1/2 


Even 


Three 





-1/2* 


1 


1/2 


Odd 


Three 





-1/2 


1 


1/2* 


Odd 



Table 3. Quark distributions in a spin-1 hadron through twist-three. 



4-3. Quark Distribution Functions for the Nucleon 

The distribution functions for a spin-i target deserve special attention because protons 
and neutrons are the principal targets of interest. In Table § the quark distribution functions 
for a nucleon target are listed through twist-three. They are classifled according to their twist 
(or light-cone projection) and their helicity. 

The distribution functions /i, gi and qt are familiar because they can be measured in 
lepton scattering. The others are less well known, but are essential to understand the nucleon 
spin substructure in deep inelastic processes. All of them are deflned by the matrix elements 
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t Wlbt LJyi/ ) 


IN diiie 


neiicity /inipiituQe 


iviecLoui enieiK 


v^'nii dii xy 


Two 


Jl 


^1 i 1 1 -|- yii 1 1 1 

-^2' 2 2' 2 


opiii civeicige 


Even 


1 WO 


9i 


^1 ill — J\^ 1 1 1 

-^2'-^2 2'-^ 2 


XlcllClLy QliicI cllCc 


Even 


Two 


Ail 


i 1 1 

^2 '-^ 2 




\JkX\1 


Three 


e 


vA-i 1 • -I 1 + ^1 1 1 1 

±2 >J-2 2' 2 


Spin average 


Odd 


Three 


hL 


^1 1 • -1 1 — ^1 1 1 1 

±2 ,-1-2 2'-^ 2 


Helicity difference 


Odd 


Three 


9t 


^^2 '-^ 2 


Helicity flip 


Even 



Table 4. Nucleon structure functions classified according to their twist and target helicities. 



of quark bilocal operators 
dX 



Ji\x 

^e^"^ l^PS\^{^)^^^^^[\n)\PS) = 2 {gi{x)p^S ■ n + gT{x)S^^ + M^g:,{x)n^S ■ n} , 

^e'^'^ (PS\^{0)ij{Xn)\PS) = 2e(x), 

^e'^- {PSmO)ia^,j,i;i\n)\PS) = 2{h,{x){S^,p, - S^^p,) / M 

+ hL{x)M{p^n^ - Pyn^)S ■ n 

+ h:i{x)M{S^^n,-S^,n^)}, (4.4) 

Some twist-four distributions (/4, (73, and h^) appear in these matrix elements. However, 
they are joined by many other quark-quark and quark-gluon distributions from which they 
cannot be separated, so there is no point in keeping track of them in this analysis. 



4.3.1. Nucleon Spin Structure at Twist-Two 

/i, gi and hi are twist-two, i.e. they scale modulo logarithms. They can be projected out of 
the general decompositions, eq. ([4. 4]), 

r d\ 

/i(^) = j ^e''%P\mm^n)\P) 
r d\ 

9i{x) = j -e''%PS\\\mil^^{^n)\PS\() 
r fix 

hi{x) = J —e'''^{PSM0)[^^,4]^,i;{\n)\PS^) (4.5) 
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To understand their physical significance — in particular, to see why a third quark 
distribution in addition to /i and gi is necessary to describe the nucleon's quark spin sub- 
structure at leading twist in the parton model — it suffices to decompose them with respect 
to a light-cone Fock space basis. If we use the helicity basis, then 

fi(x) = -(P\RUxp)R(xp)+L^{xp)L{xp)\P) 

X \ ' 

gi{x) = — {Pe^\R\xp)Rixp) — L\xp)Lixp)\Pez 

X ^ ' 

hi{x) = -Re(^Pei\L\xp)R{xp)\Pei) , (4.6) 



in analogy with eq. (|3.49|) , where we have integrated out the dependence on transverse 
momentum. Here 63 and Ci are unit vectors parallel and transverse, respectively, to the target 
nucleon's three-momentum. Clearly fi and gi can be interpreted in a probabilistic way: /i 
measures quarks independent of their helicity and gi measures the helicity asymmetry. But 
hi does not appear to have a probabilistic interpretation, instead it mixes right and left 
handed quarks. 

If instead we use a transversity basis, diagonalizing 7^, we find, 
fi{x) = —{P\a\xp)a{xp) + (3\xp)l3{xp)\P) , 

X 

2 

gi{x) = -Re{Pes\a^{xp)P{xp)\Pe3), 

X 

hAx) = —(Pei\aUxp)a{xp)—l3Uxp)(3{xp)\Pei). (4.7) 

X 

Clearly hi can be interpreted as the probability to find a quark with spin polarized along the 
transverse spin of a polarized nucleon minus the probability to find it polarized oppositely. 
fi still has the same interpretation, while now gi lacks a clear probabilistic interpretation. Of 
course the structure here is merely that of a 2 (g) 2 - spin density matrix, with the assignments 
/i ^ 1, gi *^ (J3, and hi ^ oi in the basis of hehcity eigenstates. The remaining element, a^-, 
is related to hi by rotation about the es-axis. In non-relativistic situations, spin and space 
operations (Euclidean boosts, eic.) commute and it is easy to show that gi = hi, so hi is a 
measure of the relativistic nature of the quarks inside the nucleon. 

The chirally odd structure functions like hi fig. (fig. |?I|a) are suppressed in DIS. The 
dominant handbag diagram fig. (fig. |21[b) as well as the various decorations which generate 
logQ^ dependences, aQCDiQ"^) corrections and higher twist corrections, examples of which 
are shown in figs. [21|c-e involve only chirally-even quark distributions because the quark cou- 
plings to the photon and gluon preserve chirality. Only the quark mass insertion, fig. (^Tj f), 
flips chirality. So up to corrections of order nig/Q, /ii(x, Q^) decouples from electron scat- 
tering. 

There is no analogous suppression of hi{x, Q^) in deep inelastic processes with hadronic 
initial states such as Drell-Yan. The argument can be read from the standard parton diagram 
for Drell-Yan, fig. (^2]) . Although chirality is conserved on each quark line separately^he two 
quarks' chiralities are unrelated. It is not surprising, then, that Ralston and Sopera found 
that hi{x,Q'^) determines the transverse-target, transverse-beam asymmetry in Drell-Yan. 
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(f) 

Fig. 21. Chirality in deep inelastic scattering: a) Chirally odd contributions to hi{x); b)-e) Chirally even 
contributions to deep inelastic scattering (plus L ^ i? for electromagnetic currents); f) Chirality flip by 
mass insertion. 
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Fig. 22. Chirality in Drell-Yan (plus L ^ R) production of lepton pairs . 

4.4- Transverse Spin in QCD 

The simple structure of eqs. ( [4.6|) and ( [4.7| ) shows that transverse spin effects and 
longitudinal spin effects are on a completely equivalent footing in perturbative QCD. On 
the other hand, hi was unknown in the early days of QCD when only deep inelastic lepton 
scattering was studied in detail. 

Not knowing about hi, many authors, beginning with Feynman,@ have attempted to 
interpret qt as the natural transverse spin distribution function. Since qt is twist-three 
and interaction dependent, this attempt led to the erroneous impression that transverse 
spin effects were inextricably associated with off-shellness, transverse momentum and/or 
quark-gluon interactions The resolution contained in the present analysis is summarized in 
Table ^ where the symmetry between transverse and longitudinal spin effects is apparent. 
Only ignorance of hi and h^ prevented the appreciation of this symmetry at a much earlier 
date. 





Longitudinal 
Spin 


Transverse 
Spin 


Twist-2 




hiix,Q') 


Twist-3 


hL{x,Q') 





Table 5. The transverse and longitudinal spin distribution functions through twist-three. 

Since experiments to measure hi are being planned, now is the time for theorists to 
make predictions. At this time, however, not much is known about either the general behavior 
of hi or its form in models. Here is a summary, presented in parallel with gi for the purpose 
of comparison. 
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Inequalities: 

\gi{x,Q^)\ <fiix,Q') 

\hi{x,Q')\ </i(x,Q2) (4.8) 

for each flavor of quark and antiquark. These follow from the positivity of parton 
probability distributions (see eqs. ( [4 .61 ) and ( ^4.7| )). Another inequality, proposed by 
SoffeiS has attracted attention recently, 

f^ix,Q') + g1ix,Q') >2\hUx,Q% (4.9) 

valid for each flavor (a) of quark and antiquark. Soffer's inequality is invalidated by 
QCD radiative corrections,0 in much the same way as the Callan-Gross relation, F2 = 
2xF2. Despite this problem, the inequality may prove to be a useful qualitative guide 
to the magnitude of /ii. A recent discussion of QCD radiative corrections to Soffer's 
inequality may be found in.0 

Physical interpretation: The structure function hi{x,Q'^) measures transversity. It is 
chirally odd and related to a bilocal generalization of the tensor operator, qa^ui'j^q. On 
the other hand gi{x,Q'^) measures helicity. It is chirally even and related to a bilocal 
generalization of the axial charge operator, 57^75^. Although hi is spin-dependent, it 
is not directly related to the quark or nucleon spin. It would be very useful to have a 
better idea of the dynamical and relativistic effects which generate differences between 
Qi and hi. 

Sum rules: If we define a "tensor charge" in analogy to the axial vector charge measured 
in /3-decay, 

2S'5q'^m ^ (P^|ga^^5^g|^jP5), (4.10) 

where A*^ is a flavor matrix and is a renormalization scale, then Sq"'{Q^) is related 
to an integral over h^^x, Q^), 

Sq^iQ^) = C dx{hl{x, Q2) _ hlix, Q2)) (4.11) 
Jo 

where /i" and h^ receive contributions from quarks and antiquarks, respectively. The 
analogous expressions for gi{x,Q^) involve axial charges, 

2s'Aq^{Q') ^ {PS\q^^^,^q\^jPS) (4.12) 
Aq^iQ') = I' dx{gl{x,Q^)+gl{x,Q^)). (4.13) 

Note the contrast: /ii(x,Q^) is not normalized to a piece of the angular momentum 
tensor, so /ii, unlike gi, cannot be interpreted as the fraction of the nucleons' spin 
found on the quarks' spin. Note the sign of the antiquark contributions: dq"" is charge- 
conjugation odd, whereas Aq"" is charge-conjugation even. gets no contribution from 
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quark- ant iquark pairs in the sea. All tensor charges Sq"" have non- vanishing anomalous 
dimensions, but none mix with gluonic operators under renormalization because they 
are chirally odd and gluon operators are even. In contrast, the flavor non-singlet axial 
charges, Ag'*, a 7^ 0, have vanishing anomalous dimensions, whereas the singlet axial 
charge Ag'' has an anomalous dimension arising from the triangle anomaly.E3 

• Evolution: It is worth re-emphasizing that hi has unusual QCD evolution properties. 
All of the local operators associated with hi have non-vanishing leading order anoma- 
lous dimensions. On the other hand, no gluon operators contribute to hi in any order, 
because hi is chiral odd. So hi is a non-singlet structure function - it evolves homoge- 
neously with Q^, but none of its moments are independent. 

• Models: hi and gi are identical in non-relativistic quark models, but differ in relativistic 
models like the bag model - see fig. (PB]). 

4-5. Twist-Three: Physics with g2{x,Q'^) 

There are several reasons to be particularly interested in the transverse spin dependent 
structure function, g2{x,Q'^) 

1. it can be measured in deep inelastic scattering; 

2. it is unique among higher twist distributions in that it dominates the cross section in 
a specific kinematic domain — at 90° all twist-two effects decouple, see eq. ( p,.22| ); 

3. it is related to interesting quark gluon matrix elements; 

4. it should obey an interesting sum rules. 
For a review of the properties of g2, see refill 

4.5.1. Operator Product Expansion 

The spin structure functions gi and g2 parameterize the antisymmetric part of the hadronic 
tensor Wf^^ as shown in eq. ( |1.15| ). Applying the methods of §2 we can relate the moments of 
the antisymmetric part of PF^i, to the matrix elements of quark operators accurate through 
twist-three. Consider the antisymmetric part of the forward Compton amplitude Tfj^^: 



The leading twist (twists-two and -three) contributions come from the series of operators 




2 [-^M!^ -^'^m] ' 
{PS\T[^,]{q)\PS) , where 




(4.14) 



( 



-2 



) 



n+l 



Cn{q'). 



(4.15) 



n=0,2,4... 
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Fig. 23. Bag model calculation for hi and gi from.i'i 
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with 

0a{Mi. = ^"V'7<t75^{/.i • • • Df^^yij - traces, * (4.16) 

Here { , } denotes symmetrization and [ , ] denotes antisymmetrization of the enclosed 
indices. The string of symmetric indices originates in the expansion of the bilocal operator 
about ^ = 0, and contraction of each index with q^^ selects the totally symmetric part. On the 
other hand, the index on the 7-matrix is not symmetrized. C„(g^) is the coefficient function 
and is unity to lowest order. Flavor factors such as squared charges are suppressed. 

Eq. ( |4.16|) can be split into a totally symmetric part and a part of lower symmetry, 

where the first term has twist-two (dimension n + 3, spin n + 1): 



6; 



and the second term, with no totally symmetric part, has twist-three (dimension n + 3, spin 
n): 



0[cr{^l]/X2.../t„} ^ _|_ ^ |0(T{/ii.../i„} 

+ 0<7{A»l.../i„} - &tl2{t^ia...fln} + •••}• (4.19) 

The proton matrix elements of these operators are 

PS\Q{^,,...,,^y\PS) = {S^P,,...P,^+S,,P^...P,^ + ...} 



n+1' 



PS\'S)[a{lJ.i]n2---t^n}\PS) — {('S'o--Pfli 5*^1 -Per ) -P^2 ••• -^Mn 

+ {S^P^^-S,,P^)P,,...P,^ + ...}^. (4.20) 
4.5.2. Wandzura - Wilczek Decomposition of g2 

Extracting the relation between these matrix elements and the moments of gi and g2 is an 
exercise in the methods of §2. Substituting these matrix elements into the definition of Tj^jy], 
writing a dispersion relation in terms of gi and g2 and equating terms in the Taylor expansion 
in 1/x in the BJL-limit, we find 

1 



J dxx''gi{x,Q'^) = \in n = 0,2,4. 





1 



/I Ti 
dxx^g2{x,Q'^) = (c/„-a„)), n = 2,4..., (4.21) 

4 n + 1 



*The terms denoted "traces" are whatever is necessary to subtract from the displayed term in order to render 
the resulting operator traceless and will be suppressed in the following. 
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Notice that the same operators which determine gi make an appearance in the moments of 
g2. It follows that g2 can be decomposed into two parts, one which is fixed by gi, and another 
— the "true twist-three" part — associated with the operator of mixed symmetry, 

g2ixQ^) = -giix, Q^) + / %r{y, Q^) + U^, Q'). (4.22) 

y 

X 

Wandzura and Wilczek proposed this decomposition in 1977.i^ They went further and sug- 
gested that g2 might be zero. From another perspective, g2 is the interesting part of g2. It's 
moments, 

1 

/ dxx^U^, Q') = —^d4Q'), (4.23) 
J 4(n + 1) 

^ ' 

are twist-three and measure quark-gluon correlations, as was pointed out by Shuryak and 
Vainshteyn.0 They showed that the QCD equations of motion can be used to trade the an- 
tisymmetry of Q[a{|J.l]|J,2...^J.„} factors of the gluon field strength G^u and the QCD coupling 
constant g. The result is: 

^ [l=0 

+ I E ^""'V^^m-^M.(^M,+iG'^/.,+J^M,+3-^M.-l7M.75^| , (4.24) 
^ 1=0 ) 

where Gq,/3 = \eai3X(jG^'^ and Sn symmetrizes the indices /xi . . 

Eq. ( |4.24| ) is quite formidable. A simpler example might help explain how manipulation 
of the equations of motion exposes the interaction dependence of higher twist operators. 
Consider the twist-three operator, 

X^, = ^D^D,^ (4.25) 

as an example. X is clearly twist-three (dimension five, spin no greater than two). One might 
be tempted to make a "parton model" for the matrix element of X, by replacing — > 9^ 
and evaluating X in a beam of coUinear quarks. That, however, would be a mistake, since 
application of the identity = ^{jfj., Ip}, and the QCD equations of motion, Ip^p = iplp = 
and [D^, D^] = gO^^y, yields 

X^, = ^-^G^xl^-i,ij. (4.26) 

So it is clear that X and its matrix elements are interaction dependent and measure a 
quark-gluon correlation in the target hadron. This is a special case of a general result that 
operators of twist > 3 can always be written in a form in which they are manifestly interaction 
dependent. In particular one has: 

\%{x,Q^) oc {PS\^gS^,^,i>G^,,j^,^\PS). (4.27) 

Model builders or lattice enthusiasts who want to predict g2 must confront such matrix 
elements. 
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4.5.3. The Burkhardt Cottingham Sum Rule 

A striking cousequeuce of the light-coue analysis of g,_ is the appaieut sum rule.S 

1 

J dxg2{x,Q^) = 0, (4.28) 



which foUows directly from eq. ( [4.21D with n = 0. If true, the sum rule requires g2 to have a 
node (other than at and 1). 

The Burkhardt-Cottingham (BC) Sum Rule looks to be a consequence of rotation 
invariance. To see this, return to eq. ( [4.4| ) and consider the 7^75 case in the laboratory frame 
(where = 0). First, set S = Mis, contract the free Lorentz index with and integrate 
overall x, leaving 

1 

J dxgi{x,Q') = (Pe3|g(0)7Sg(0)|Q2|Pe3). (4.29) 
-1 

Next repeat the process with S = Mii and /i = 1, with the result, 

1 

J dxgT{x,Q') = (Pei|g(0)7S5g(0)|Q2|Pei). (4.30) 
-1 

The right hand sides of these two equations are equal in the rest frame by rotation invariance, 
so 

J'^ dxgrix, Q2) = J'^ dxg,{x, Q^), (4.31) 

whence ^ 

J ^dxg2{x,Q^) = 0, (4.32) 

apparently a consequence of rotation invariance. 

The subtlety in this derivation is that the integral in eq. ( [4.32 ) goes from — 1 to 1 
including x = 0. As we have defined it, g2{x,Q'^) is the limit of a function of and u 
and therefore might contain a distribution ((5-functions, etc.) at x = 0. Suppose g2 has a 
5-function contribution at x = 0, 

92{x, Q') = gt''"^'"%x, Q') + c5{x). (4.33) 

Then since experimenters cannot reach x = 0, the BC sum rule reads 



/ dxg. 
Jo 



1 

— ( 
2 



observable^^^g2) _ (4.34) 



which is useless. 

This pathology — a 5-function at x = — is not as arbitrary as it looks. Instead it 
is an example of a disease known as a " J = fixed pole with non-polynomial residue". First 
studied in Regge theory, a 6{x) in g2{x, Q"^) corresponds to a real constant term in a spin 
flip Compton amplitude which persists to high energy. There is no fundamental reason to 
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exclude such a constant. On the other hand the sum rule is known to be satisfied in QCD 
perturbation theory through order 0{g'^). The sum rule has been studied recently by several 
groups who find no evidence for a 5{x) in perturbative QCD.&S So at least provisionally, we 
must regard this as a reliable sum rule. At least one other sum rule of interest experimentally, 
the Gerasimov, Drell, Hearn Sum Rule for spin dependent Compton scattering has the same 
potential pathology. For further discussion of the BC sum rule see ref..li^ 



4.5.4. The Evolution of g2 

The following discussion can be regarded as a "theoretical interlude" and may be omitted 
by the casual reader. 

The evolution of quark distribution functions is an unavoidable complication in per- 
turbative QCD. Data are inevitably taken at different Q^, making it difficult to evaluate sum 
rules (which some definite Q"^) without some information about Q^-evolution. 

At leading twist the subject is well understood and we have ignored it. The evolution of g2 
is more complicated. Eq. ([4.4]) provides as deceptively simple operator representation of qt, 



^e'''{PS^\m^±l^H>^n)\Q.\PS,_) (4.35) 

It looks as though grp is determined by a single operator which evolves homogeneously. 
However, the operator in eq. ( |4.35| ) is equivalent to the series of quark gluon operators given 
in eq. ( [4.24 ). Worse still, all these operators mix under renormalization. The number of 



operators that mix grows linearly with n. Back in the 1970's, Ahmed and Ross calculated 
the anomalous dimension matrix element for the evolution of the transverse quark axial 
current operator in eq. ( [4.35|) .il In 1982 Shuryak and Vainshteyn pointed out that this was 
but one element of an x anomalous dimension matrix M Later Ratcliffe, Lipatov et al. 
and others calculated the full anomalous dimension matrix.&i3 

Underlying this matrix renormalization group evolution is a somewhat simpler parton 
picture. The fundamental objects for studying the evolution of g2 are two-variable parton 
distribution functions defined by double light-cone fourier transforms.! 



/— ^e^'^'+^'^(^^-^')(P5|g(0)iZ^"(^n)7^g(An)|Q.|P5) = 2ie"^'^^npS,p,G{x,y,Q^) 

+ ... 

/ ^^e^'''"'^(^"'^(^^l^(0)^^"(/^^)'^75g(An)|Q2|P^) = 2ie-P"-npS„PrG{x,y,Q^) 



+ ... 

(4.36) 

where . . . denote other tensor structures of twist greater than three. The functions G{x^ y, Q^) 
and G{x, y, Q"^) are generalizations of parton distributions describing the amplitude to find 
quarks and gluons with momentum fractions x, y and x — y m the target nucleon. They can 
be represented diagrammatically as in fig. (2^). These "higher twist distribution functions" 



share many properties with leading twist. Only good light-cone components and coUinear 
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Fig. 24. Quark/gluon diagram contributing to the two-variable parton distribution functions G and G. 

momenta appear - transverse momentum has been eliminated in favor of the interactions 
which generate it. The variables x, and x — y take on only physical values, —1 < (x, y,x — 
y) < 1, corresponding to emission or absorption of quarks, antiquarks or gluons.0 

At any value of Q^, grix, Q"^) can be projected out by integrating out the variable 

9t{x, Q2) = 1 j'^ dy{G{x, y, Q') - G{y, x, Q^) + g(x, y, Q^) _ G{y, x, Q^) } (4.37) 

The reason for introducing G and G is that they obey natural generalizations of the Gribov- 
Lipatov-Altarelli-Parisi (GLAP) evolution equations. Schematically, 
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(4.38) 



for some "splitting function" , V . In order to evolve from one to another it is necessary 
to know G and G as a function of both x and y. However, measurement of gT does not 
supply that information, is, in a sense, a "compressed" structure function. Information 
essential to evolution has been integrated out of it. The associated "open" distributions, G 
and G evolve simply. gT does not. 

The situation is not entirely hopeless, however. It aopears that the evolution of g2 
may simplify in the A^^c oo limit, especially at large x.^ Ali, Braun and Hiller showed 
that in this limit g2 obeys a standard GLAP equation, albeit with non-standard anomalous 
dimensions. As measurements of g2 improve, theorists will be forced to expend more effort 
to clarify the nature of evolution of gT- 



4.5.5. The Shape of g2 

With the first measurements of g2 now available]^ and more accurate measurements expected 
soon, it seems timely to review model predictions. First, here is an agenda of progressively 
more detailed questions for experimenters as they relate to theory 

• Is g2 zero? 

The answer is in on this one: the new SLAG data show (72 7^ 0. 
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• Is Jq dxg2{x, Q^) zero as required by the BC sum rule? 

The SLAG data are consistent with the BC sum rule, albeit with large errors. 

• If Jo dxg2{x,Q'^) = 0, then g2 must have at least one non-trivial zero (besides x = 1 
and perhaps x = 0). Is there just one such node? 

The SLAC data are consistent with one node, again within large errors. 

• If there is only one non-trivial node, what is the sign of 



With only one non-trivial node, the sign of M2 determines the gross structure of g2- If 
M2 is positive g2 is negative at small x and positive at large x, and vica versa if M2 is 
negative. The Wandzura-Wilczek (WW) contribution to (72 gives < 0. 

The SLAC data favor M2 < 0. 

• Is there a signal for ^2 or does the WW contribution account for all of (72? 

The SLAC data agree rather well with the WW prediction, however the accuracy is 
such that a fairly significant g2 term would not yet be detectable. Some observers claim 
to be able to see a deviation from WW with M2 > 0. 

• As more accurate data become available it will be possible to subtract away the WW 
contribution to reveal g2- 

g2 has been studied in quark models and using QCD Sum Rules. Perhaps the most 
thorough analysis, including QCD evolution, has been performed by Strattman.ii He has 
taken several versions of the MIT bag model, calculated g2 and then evolved the resulting 
distribution using the method of Ali, Braun and Hiller to experimentally interesting Q^. His 
estimates for the proton are shown in comparison with the SLAC data in fig. (p5|). 

His estimates of ^2 are small compared to g^^ and cannot be excluded by the existing 
data. The neutron g2 is very small in quark and bag models for the same reason that the 
neutron's gi is so small: the correlations of charge and spin in the SU (6) symmetric neutron 
tend to cancel for the neutron. The second moment of g2 has also been estimated using 
QCD sum rules.il Surprisingly ref.@ finds M2 [(72'^°*°'^] consistent with zero and M2 [(72^"*'^°'^] 
negative and significantly different from zero. Existing data cannot rule out this behavior — 
we shall have to wait for the HERMES facility at HERA to provide more accurate data. ^2 
depends on quark gluon correlations within the nucleon, which are not likely to be perfectly 
described in such simple models, so these predictions should be more as a guide than a 
prediction of expected behavior. 

5. The Drell-Yan Process 

Although deep inelastic scattering has been the source of much insight into nucleon 
structure, it has many limitations: polarized targets can only be probed with electromagnetic 
currents (neutrino scattering from polarized targets being impractical); gluon distributions 




(4.39) 
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Fig. 25. Data on (72 from the SLAC E143 collaboration compared with the WW contribution and the 
bag model estimates of Strattmann. 
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do not couple directly, but instead must be inferred from careful study of the evolution of 
quark distributions; chiral odd distributions like hi decouple. None of these limitations af- 
flict deep inelastic processes with hadron initial states. These include not only the original 
Drell-Yan process, pp ^ i£ + X via one photon annihilation, but also generalizations to 
annihilation via and Z^, and parton-parton scattering resulting in jet production or 
production of hadrons or photons at high transverse momentum. QCD predictions for all of 
these processes are obtained by combining quark/gluon distribution and fragmentation func- 
tions with hard scattering amplitudes calculated perturbatively. This formalism is treated 
in standard texts — our object is to explain how the spin, twist and chirality classification 
developed in previous sections can be applied to Drell-Yan processes. In this section we will 
treat only the "classic" Drell-Yan process, pp —>■ 'j* + X ii + X . The generalizations of our 
spin/twist /chirality analysis to other processes is fairly straightforward and yields interesting 
predictions for a variety of processes. The original treatment of polarized Drell-Yan at this 
level was made by Ralston and Soper.il 

When last we considered Drell-Yan in §2.4, we noted that the dominance of the leading 
light-cone singularity could be overwhelmed by rapid phase oscillations in the matrix element. 
Fig. (|r^) shows a contribution to Drell-Yan which is obviously proportional to the product of 
quark distribution functions each of which has the form of eq. ( p.32|) . Each oscillates rapidly 
along the light-cone oc e^°'^'^, where P is either of the two external hadron momenta. In this 
section we will use the formalism of the previous two sections to compute fig. ([TTD , to show 
that it generates large (scaling) contributions to Drell-Yan, and to classify them with respect 
to spin and twist. 

5.1. Operator Analysis 

We begin with the Drell-Yan tensor, Wfj_^, from eq. (|2.29| ), with momenta and spin 
more carefully labeled (and IN labels suppressed): 

W,. = ^J d^ie^'-HPASAPBSBVJMUOlPASAPBSB). (5.1) 

The dominant contribution for the Drell-Yan process is shown in fig. (|26|). From the diagram 
it seems that W^j^u reduces to a product of light-cone quark correlation functions, as indicated 
by the markings on the figure. The essential features are 1) that nothing propagates between 
the two currents, so there is no singularity as — 0; and 2) each current has a quark line 
landing on each hadron. The latter suggests that the diagram can be factored into products 
of quark-hadron amplitudes by making a Fierz transformation in order to couple the spin, 
color and flavor indices on the quarks in a more appropriate order. First we write out the 
currents in terms of quark fields, limiting ourselves to terms symmetric in /i z/ (which 
survive contraction with the lepton tensor, L^^), 

j{M(o)J-}(0 = -^fc(0)V'z(OV^.(OV^.(0)(l7^'^)fc,(l7''^).z. (5.2) 

Note the — sign from anticommuting quark fields. We wish to recouple indices so quarks 
acting in the same hadron are coupled to one another — {kj){il) {kl){ij). The color Fierz 
transformation is simple, 

l,,l,, = ^(l«li, + 2^^), (5.3) 
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where 1 is the 3x3 unit matrix in color space. [This resuh is easily derived by multiplying 
both sides by the matrices Ikj and \kj, in turn, and using multiplication properties of the 
A's.] The second term in eq. (|5.3|) vanishes when matrix elements are taken in color singlet 
hadron states. The Fierz transformation for the Dirac matrices is more complicated but 
better known, 



(5,4) 



where the omitted terms are traces of the terms shown explicitly. One is left with bilocal 
light-cone correlation functions of the form: vector x vector, axial vector x axial vector and 
tensor X tensor. The flavor Fierz transformation is straightforward and is left to the reader. 

After reorganizing the product of currents to group fields that act in the same hadron 
together, the matrix element in eq. ( |5.1| ) factors into the product of quark distribution 
functions. However, the coordinate interval is not constrained to be lightlike. In ref.i 
each bilocal matrix element is expanded about the tangent to the light-cone defined by the 
large momentum or Pg, and it is shown that the light-cone contributions dominate. The 
matrix element in hadron A is a non-trivial function of Pa ■ ^ along the surface = and 
Pb ■ ^ = 0, while the matrix element in hadron i? is a non-trivial function of Pb ■ ^ along 
the corresponding tangent to the light-cone. We refer to ref.i for further details. Accepting 
this, Wfj_i, reduces to the product of quark distribution functions, one for hadron A, another 



for hadron B. The VxV,AxA, TxT structure of eq. (|5.4| ) is reflected in the structure 
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of the resulting product of distribution functions. The V x V part is described by: 



= ^{2nY6\Q^) E elf^{x)ff{y) Wb + - ^Va ■ Pb) , (5.5) 

"J a 

where 

Pi^ = p'a + —p''b (5.6) 
s 

= p'b + ^pI, (5.7) 

with p\ = Pb = and ■ pb = s. This result is valid up to corrections of order or 
1/s. There are no order 1/Q or l/^/s (twist-three) corrections — another example of the 
general result that only even twists appear in spin average deep inelastic phenomena. The 
5(Q5_) in sq. (|5.5|) reflects the fact that to leading twist and leading order in QCD radiative 
corrections, it is as though all partons move parallel to the parent hadron. Our calculation 
can only be used to study observables that integrate over Q^. Eq. (|5.5| ) is the original result 
of Drell and Yan: quarks that annihilate to form a virtual photon of squared-mass and 
three-momentum must have xy = Q'^/s and x — y = 2Q^/y/s. 
In the same spirit the A x A contribution is 

= -]:i2nmQ^) I J: elgl{x)gl{y f^- ^^^^ "f ^ (p>^ + - g^^p^ ■ ps) , 
3 \a [Pa- PbY 

+ Y.^l9l{x)g^{yf-^{p^^S''B^+f^S%^) (5.8) 

a PA-PB 
+ Eea^7T(^M(l/)^^(pS^AX+PB^^j|, 

a Pa- Pb J 

The first line is twist-two and contributes only when both initial hadrons are longitudinally 
polarized. Not surprisingly, this contribution measures gi®gi- The latter two lines are twist- 
three, suppressed by the factor S'_l, and contribute only when one hadron has longitudinal 
polarization and the other transverse. It is worthwhile relating the spin, twist and chirality 
structure of this result to the classification scheme developed in §'s 3 and 4. Products of the 
form gi ® gi and gi ® gr conserve quark chirality and contribute at orders of l/Q"^ which 
reflect the twist assignments we made in §4. Axial vector bilocals can only generate gi and 
gTi so we can be confident that we have not missed other contributions. In fact, this result 
could have been written down a priori, up to coefficients of order unity simply by carefully 
considering the selections rules and twist assignments developed earlier in these notes. What 
is not at all obvious, however, is that other twist-three quark/gluon operators, of the form 
discussed in the previous section, and not directly related to gr, do not arise when gluonic 



corrections to the diagram of fig. ^ are computed. We take up this question below. 



Finally, the T xT contribution takes the form, 

= -l{27r)'5\Qj[SA±-SB±ip'APB+PAP'B-9''''PA-PB) 
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+ iPA ■ pb) (^^x^^x + s'-^^s'sx) 1 E ^iKi^my) 



- T.elKi^)hliy)^^-^ip'AS'l^+PAS'A^) (5.9) 

a Pb - Pa 

- T.^lhli^)K{y)^^^{p'BSB^+PBS'B^)]- 

a PA-PB J 

Here the first two lines are twist-two — they scale (modulo logarithms) in the deep inelastic 
limit. They contribute only when both initial hadrons are transversely polarized and provide 
a leading twist probe of transversity distributions. The second two lines are twist-three 
and contribute when one hadron is transversely polarized and the other longitudinally. Once 
again the classification scheme of §'s 3 and 4 is illustrated. Note, for example, that transverse- 
longitudinal polarization receives contributions from both gi ® Qt and hi^h^. 

Unfortunately things are not quite as simple as they have been presented so far. Two 
questions of gauge invariance arise — one straightforward and the other rather subtle: 

• The factorization of W^^ into products of bifocal operators acting in different hadrons 
is not color gauge invariant. Evidence of this is the absence of the Wilson link, 

v{e^Vi f^dCA.iC)^ (5.10) 

between the quark fields at and ^. The problem is that we have not included those 
diagrams that represent each quark propagating in the color field of the remnants of 
the other hadron. The same problem was discussed in §3 for the case of deep inelastic 
scattering, where it is easier to handle because the operator product expansion will 
preserve gauge invariance if treated with sufficient care. In this case we must find and 
analyze the diagrams which restore color gauge invariance through twist-three and 
be certain that they do not generate any contributions to W^iv beyond rendering the 
identification of Drell-Yan matrix elements with quark distributions functions gauge 
invariant. An example of the type of diagram that does the trick is shown in fig. (|27|) . 
The interested reader is referred to ref.i for details. 



The results we have just quoted (eqs. (|5.9|) , (|5.5| ), and ( |5.10|) ) violate electromagnetic 



gauge invariance at twist-three. It is easy to see that the twist-three terms do not satisfy 
WiJLv(f = 0. Once again the answer lies in diagrams like fig. (0). In this case, when 
the Wilson link which appears in the bilocal operator is expanded to first order away 
from the light-cone, a set of contributions with explicit transverse gluon fields arise. 
Using the equations of motion these can be related back to the product of good and 
had quark fields which define the twist-three distributions qt and h^. With sufficient 
care one finds exactly the terms necessary to restore electromagnetic gauge invariance. 
Once again a fuller discussion can be found in ref..i 
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Fig. 27. New class of diagrams, needed to restore color and electromagnetic gauge invariance. 



5.2. Polarized Drell-Yan: a Brief Summary 

Since the equations and the analysis in this section have become rather comphcated, 
it is useful to extract the simple predictions for spin asymmetries to provide a summary and 
to reference the workers who originally derived these results. 

In the polarized Drell-Yan process, three cases appear: longitudinal-longitudinal {LL) , 
transverse-transverse (TT) and longitudinal-transverse {LT), as illustrated in fig. (^). If 
instead of virtual photon production, we had considered Drell-Yan production of or 1^^, 
then a longitudinal, leading twist, parity violating single spin asymmetry proportional to 
Qi ® fi appears. There is no analogous transverse asymmetry because the product hi ® fi 
cannot conserve ouark chirality. The longitudinal asymmetry, which was first studied by 
Close and Sivers,E3 is given by: 

Ralston and Soper first discovered transversity, defined a twist-two distribution and expressed 
the transverse asymmetry as:0 

_ sm'e cos2cPY.aelK{x)h\{y) 

The angles are defined in the lepton center of mass frame. The longitudinal-transverse asym- 
metry has been investigated by Jaffe and JiH and can be written as: 

^ 2sin2gcos0 M EaeMix)ygUy) - xhl{x)hliy)) 

l + cos2^ Eae^/f(x)/f(y) • ^ ■ ' 

Clearly, it is a twist-three observable, which in principle allows for a measurement of hi- 

6. Annihilation and Quark Fragmentation 
Functions 

As a final application we give a brief introduction to the classification and uses of 
the spin dependent fragmentation functions which determine the distribution of final state 
hadrons in deep inelastic processes. There are strong reasons to want to develop a better 
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Fig. 28. Various polarization configurations, studied in the Drell-Yan process. The small arrows denote 
the spin projections of the two nucleons relative to the direction of the nucleon beams (long arrows). For 
one of the nucleons the difference of the up and down projections is taken. One distinguishes longitudinal 
(top), transverse (middle) and longitudinal-transverse (bottom) scattering. 

understanding of hadron fragmentation processes. In §1 we mentioned the possibility of 
studying the spin structure of unstable hadrons like the A, p and D* . Another reason is 
that parity violating processes like qq ^ hadrons provide probes of spin structure 

unavailable in deep inelastic scattering, where the analogous experiment would be neutrino 
scattering from a polarized target. A final reason is that the selection of a particular hadronic 
fragment can serve as a filter for an interesting quark or gluon distribution function. The 
material in this section is based primarily on refs.0il andil where more details can be found. 

6.1. Single Particle Inclusive Annihilation 

The simplest quark fragmentation function is represented diagrammatically in fig. (|29|) . 
More complicated fragmentation processes, such as coherent fragmentation of several quarks 
and gluons, do not contribute until order beyond our interest. First we consider the 

helicity classification in analogy to §4.1-4.3. In the figure a quark of momentum k and helicity 
A fragments into a hadron of momentum P and helicity A plus an unobserved final state 
X. The process then repeats in reverse as the unobserved system, X, plus the hadron of 
momentum P and helicity A' reconstitute the quark of momentum k and helicity A'. The 
scattering k + P^k + P is forward, i.e. collinear. For definiteness, we take the momentum 
of the quark-hadron system to be aligned along the es-axis. Then helicity is conserved as 
a consequence of angular momentum conservation about this axis: A — A = A' — A'. As in 
deep inelastic scattering, the initial and final hadron helicities A and A' need not be equal 
because the hadron need not have been in a helicity eigenstate; likewise for the quark. As in 
scattering, the quark lines may correspond either to good or bad light-cone components. 

Some of the results of refs. §andi are as follows. Quark fragmentation functions of 
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Fig. 29. Quark fragmentation function in a helicity basis. 



the form shown in fig. (^) and the equivalent gluon fragmentation functions (without fur- 
ther active parton hues) are sufficient to characterize hadron production in hard processes, 
provided: i) one studies leading twist ((9(1/Q°)) in any hard process, or ii) one studies an 
effect in deep inelastic lepton scattering at the lowest twist at which it arises, and one ignores 
QCD radiative corrections.^ Each appearance of a bad component of the quark field costs 
one power of \fQ^ in the deep inelastic limit [i.e. it increases the twist by unity). As in 
scattering, for produced hadrons of spin-1/2, helicity differences are observed in longitudinal 
spin asymmetries; helicity flip is observed in transverse spin asymmetries. Since perturbative 
QCD cannot flip quark chirality (except through quark mass insertions which we assume to 
be negligible for light quarks), chirally-odd quark distribution and fragmentation functions 
must occur in pairs. 

Fragmentation functions can be labeled uniquely by specifying the helicity of quarks 
and hadrons and the light-cone projection of the quarks in direct analogy to the classifica- 
tion of distribution functions. We denote fragmentation functions in a helicity basis by A. 
Parity invariance of QCD requires: ^aa,a'A' = -A_a-a -a'-A'- Time reversal invariance, which 
further reduces the number of independent quark distribution functions does not generate 
relationships among the {^1} because it changes the o'ui-state (PX)out in fig- ( P^ to an 
m-state. As in the scattering case, we denote the appearance of bad light-cone components 
by an asterix on the appropriate helicity index. 

As a simple example, consider production of a scalar meson like the pion. Through order 
l/\/Q^ there are three independent fragmentation functions: in, -^in I'n! and ^i*n in- 

2 '2 2 '2 2 "'2 

The first is twist-two and scales in the oo limit, the latter two are twist-three and are 

suppressed by l/\/Q^ in the ^ oo limit. The first function, ^iq^iq, is proportional to the 
traditional fragmentation function D{z,Q'^). It has the same twist, light-cone, helicity and 
chirality structure as /i(x, Q^), so to avoid an explosion of notation we denote it by fi{z, Q"^) 
[We will follow the same convention for other fragmentation functions.]: 

A(z,Q2)oc^o,io (6-1) 

If we were studying quark distribution functions, the other two helicity amplitudes would be 
equal by time-reversal invariance. Here, there are two independent fragmentation functions. 

61(2;, oc iio,i*o + -4i*0,iO 
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(6.2) 

We have found that the hehcity classification of fragmentation functions is identical to that of 
distribution functions at leading twist. At twist-three, however, there are more fragmentation 
functions due to the absence of time reversal constraints. 

The application to spin-1/2 is analogous to the classification of spin-1/2 distribution 
functions given in Table § of §4.2.2 except that each twist-three fragmentation function 
comes in two forms, one even and the other odd under time reversal. We suspect that the 
T-even functions are more important than the T-odd since the latter vanish in the absence 
of final state interactions. So fragmentation at twist-two requires fi{z), gi{z) and hi{z) with 
helicity, transversity and chirality properties identical to the analogous distribution function. 
The only twist-three fragmentation function likely to be of much interest is the spin-averaged, 
chiral-odd, time reversal even function ei{z). 

In order to relate particular deep inelastic processes to quark distribution and frag- 
mentation functions and to study them in models of non-perturbative QCD, it is necessary 
to have operator representations for them. This formalism is developed in refs..0il Here we 
display the results for the three leading twist fragmentation functions for a spin-1/2 hadron, 
and for the twist-three spin average function, ei. Generalizations to spin-1 can be found in 
ref..0 The generic expression for a fragmentation function takes the form, 

f(^) = Tr|r„^X: / ^e-^'^^ {0\MO)\PX) {PX\MXn)\0)Y (6.3) 

where Ta/B stands for an arbitrary Dirac matrix. This result holds under the condition that 
the diagram of fig. ( pOj) dominates. To obtain /i, gi and hi one chooses F = 7:^75 and 
aln^i'')^ respectively, 

dX 
— f 
27r 

d\ 



M 



9i{z) = ^-j:J^e-^'/^0\4^,m\PSX){PSX\^iXn)\0 

= ^E/^e-^'^'(0kX^75^(0)|P5±x)(PSxX|^(An)|0 



(6.4) 



At twist-three the equations of motions can be used to express the structure functions in 
terms of independent degrees of freedom, quantized at = 0. For example one can ob- 
tain two expressions for the chiral odd, T-even spin independent fragmentation function 
ei{z). One involving ipip and another where bad light-cone components have been traded for 
transverse derivatives and gluon degrees of freedom. 



Meiiz) = lE/^e-^V.(0|^(0)|PX)(PX|^(An)|0) 
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- {Q\ij+mPX){PX\i;+{\n)ip^{\n)m}. 

(6.5) 

With these ingredients we are now prepared to explore a few apphcations of spin dependent 
fragmentation functions. 

6.2. Polarized q ^ A Fragmentation Functions from e~^e~ ^ A + X 

This section is based on work done with M. Burkardt.l2i In the symmetric quark model, 
the A-baryon has a rather simple spin-flavor wavefunction. All its spin is carried by the s- 
quark, while the wrf-pair is coupled to = 0, / = 0: Au^ = Ad^ = and As^ = 1. While 
the quark model identifies the A-spin with the spin of the s-quark, the data on the quark 
spin structure of the nucleon suggests that the actual situation is more complex. If we take 
the latest SMC/SLAC data on the fl dxg^{x), combine it with /3-decay data and assume 
exact S'f/(3) flavor symmetry for baryon axial charges, we find 

Au^ = Ad^ = ^(S-D) = -0.23 ±0.06 



As^ = ^(S + 2D) = +0.58 ±0.07 



(6.6) 



It would be exciting to test the SU (3) flavor assumption by observing deep inelastic scattering 
from a A target. Unfortunately we have to settle for observing the A as a fragment in 
annihilation processes. The parity violating, self analyzing decay of the final state A makes 
it particularly easy to study its polarization in fragmentation processes. Measurement of the 
helicity asymmetries for semi-inclusive production of A's in e~^e~ annihilation near the 
resonance allows a complete determination of the spin-dependent fragmentation functions 
for the different quark flavors into the A. In the event that these could be measured it would 
be very interesting to compare the spin fractions measured in fragmentation, As^, Ad^, and 
Ait^, with the predictions, eq. ( |6.6D in order to get a better understanding of the role of spin 
in the fragmentation process. 

We are concerned here with twist-two helicity asymmetries, described by the fragmen- 
tation function gi{z, Q"^). For simplicity we adopt a more conventional parton-model notation 
where we define d^j^^^^ to be the probability that a left handed A fragments into a left handed 
quark, etc. Then the unpolarized differential cross section for e~e~^ A ±X is obtained by 
summing over the cross sections for e+e~ — > qq, weighted with the probability dg{z, Q^) that 
a quark with momentum -P fragments into a A with momentum P. As usual, we suppress 
the dependence generated by radiative corrections in QCD, 



There is a single polarized fragmentation function for each flavor of quark or antiquark. 



Aq{z) = 
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= 4S(^)-4S(^), (6-8) 

and furthermore isospin invariance requires that Au(z) = Ad(z) and Au{z) = Ad{z), so the 
number of independent fragmentation functions is reduces to four — e.g. Au, Au, As and 
At. 

As an example we quote the prediction for A production in e~e~^ annihilation via 
photons. In this case one has to start from polarized e~ (or e^) in order to fix the polarization 
of the quarks. One finds for the helicity asymmetric cross- sect ion, 



^2^e-{L)e+^A{L)X ^e- {L)e+ ^k{R)X ^2 

■ COS 9 



^ ^Au{z) + Au{z)) + ^ (As{z) + As{z 



9 V ^ ^ ^ V 9 



dVtdz dVLdz 2s 

(6.9) 

where L, R denotes the helicity of the e and the A. So annihilation into a single photon 
allows one combination of the four independent fragmentation functions to be measured. As 
described in ref.,il annihilation at the and just off the Z° peak where 7 — Z interference 
is largest will allow independent detection of all four quark and antiquark fragmentation 
functions. For further discussion of this process, backgrounds and experimental possibilities 
see ref.S and recent papers by the Aleph and Delphi collaborations. 

6.3. Observing hilxjQ"^) in Electroproduction 

This section is based on work done with X. Chirally odd quark distributions are 

difficult to measure because they are suppressed in totally-inclusive deep inelastic scattering. 
So far, the only practical way to determine hi{x, Q"^) we have discussed has been Drell-Yan 
production with transversely polarized target and beam. 

As an application of the fragmentation function formalism - one of many - we show 
how a chirally odd fragmentation function can be exploited to enable a measurement of 
hi{x, Q^) to be obtained in polarized electroproduction of pions from a transversely polarized 
nucleon. This is an experiment which could be performed at several existing facilities. Related 
suggestions involvingsemi-inclusive production of A-hyperons and of two pions have been 
discussed previously.OE3 The proposal outlined here is simpler since it involves only one 
particle in the final state and does not require measurement of that particle's spin. The price 
we pay for this simplicity is suppression by a power of y/Cp. 

We consider pion production in the current fragmentation region of deep-inelastic scat- 
tering with longitudinally polarized leptons on a polarized nucleon target. The simplest 
diagram for the process is shown in fig. (^OD, where a quark struck by the virtual photon 
fragments into an observed pion plus other unobserved hadrons. The cross section is propor- 
tional to a trace and integral over the quark loop which contains the quark distribution func- 
tion and fragmentation function. Due to chirality conservation at the hard (photon) vertex, 
the trace picks up only the products of the terms in which the distribution and fragmenta- 
tion functions have the same chirality. When the nucleon is longitudinally polarized (with 
respect to the virtual-photon momentum), the twist-two, chirally even distribution gi{x) can 
couple with the twist-two chirally even fragmentation function fi{z), producing a leading 
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contribution 0{1/Q^) to the cross section. On the other hand, in the case of a transversely 
polarized nucleon, there is no leading-order contribution. At the next order, the nucleon's 
transversity distribution hi{x) can combine with the twist-three chirally odd fragmentation 
function ei{z), and similarly grix) can combine with the chirally even transverse-spin distri- 
bution fi{z). Both couplings produce l/y/Q^ contributions to the cross section. It is simple 
to see, however, that fig. (0) alone does not produce an electromagnetically-gauge-invariant 
result. This is a typical example of the need to consider multi-quark/gluon processes be- 
yond twist-two. In the present case (twist-three), however, the contributions from coherent 
scattering can be expressed, with novel use of QCD equations of motion, in terms of the 
distributions and fragmentation functions defined from quark bilinears. This is a specific ex- 
ample of another rule quoted in §6.2: "Parton diagrams (without further active parton lines) 
are sufficient to characterize hadron production in hard processes, provided: one studies an 
effect in deep inelastic lepton scattering at the lowest twist at which it arises, and one ignores 
QCD radiative corrections." The combined result is gauge invariant, as can be seen from the 
resulting nucleon tensor, 

W^-^ = -ie^^-^^^^KS ■ n)ppGi{x, z) + S^fiGrix, z)\ (6.10) 

The two structure functions in W^^ are related to parton distributions and fragmentation 
functions, 



Gt(x,z) 



(6.11) 



where the summation over a includes quarks and antiquarks of all flavors. 
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To isolate the spin-dependent part of the deep-inelastic cross section we take the dif- 
ference of cross sections with left-handed and right-handed leptons, we use 



dE'dVL EM 



N 



It is convenient to express the cross section in terms of scaling variables in a frame where 
lepton beam defines the es-axis and the 61 — 63 plane contains the nucleon polarization vector, 
which has a polar angle a. In this system, the scattered lepton has polar angles (6*, 0) and 
therefore the momentum transfer q has angles (6', tt — 0). Then we obtain an expression for 
the semi-inclusive process quite similar to that for the total inclusive scattering defined by 
eq. O). 



dx dy dz dcj) Q'^ 



cos «(! — ^)G'i(x, z) 



+ cos0sina-^ (k — 1)(1 — y) ( Gt{x, z) — Gi{x, z){l — ^) 



(6.13) 



where k = 1 + 4x^M^/Q^ in the second term signals the suppression by a factor of l/\/Q^ 
associated with the structure function Gt- The existence of Gi in the same term is due to a 
small longitudinal polarization of the nucleon relative to q when its spin is perpendicular to 
the lepton beam. 

Eq. (|6l3D is the main result of this section. As a check, we multiply by z, integrate 



over it and sum over all hadron species. Using the well-known momentum sum rule, 

Y: Jdzzf^{z) = l, (6.14) 

hadrons 

and the sum rule, 

J2 J dze1{z) = (6.15) 

hadrons 

which is related to the fact that the chiral condensate vanishes in the perturbative QCD 
vacuum, we get the result for total inclusive scattering, given in eq. ( |1.22| ) The similarity 
between the inclusive and semi-inclusive cross sections suggests that they can be extracted 
conveniently from the same experiment. 

The aim of this example was to show that an unfamiliar fragmentation function (ei) 
could be employed to obtain a measurement of an interesting, if unfamiliar, distribution 
function (hi). It is apparent from eq. (|6.13|) that we have been only partially successful: 
although the hf distribution for each quark flavor appears in eq. ( |6.13| ), the sum over flavors 
couples it to the unknown flavor dependence of e". Perhaps flavor tagging can be used at 
large-z to identify the contributions of individual quark flavors. For x in the valence region 
(where one can ignore antiquarks in the nucleon), and z — > 1, the dominant fragmentation, 
u —>■ TT^, d ^ ^-.-^ K~ 1 '^^y allow one to trigger on the contributions of m, d and s 
quarks separately.L3 One might be concerned that the unknown fragmentation function, Ci, 
might not respect the dominant fragmentation selection rules, which have only been tested 
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for the spin-averaged, twist-two fragmentation function, /i. However, the coherent gluon 
interactions which distinguish the twist-three Ci from fi are flavor independent and should 
not alter the selection rules. Although this may be a difficult path to measuring hi, so are 
all the others. This one owes its existence to our improved understanding of the spin and 
twist dependence of quark fragmentation functions, including the spin-average twist-three 
fragmentation function ei. 
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